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Chapter 1 
Introduction 



In general we do not like to wait. But reduction of the waiting time usually requires extra 
investments. To decide whether or not to invest, it is important to know the effect of 
the investment on the waiting time. So we need models and techniques to analyse such 
situations. 

In this course we treat a number of elementary queueing models. Attention is paid 
to methods for the analysis of these models, and also to applications of queueing models. 
Important application areas of queueing models are production systems, transportation and 
stocking systems, communication systems and information processing systems. Queueing 
models are particularly useful for the design of these system in terms of layout, capacities 
and control. 

In these lectures our attention is restricted to models with one queue. Situations with 
multiple queues are treated in the course "Networks of queues." More advanced techniques 
for the exact, approximative and numerical analysis of queueing models are the subject of 
the course "Algorithmic methods in queueing theory." 

The organization is as follows. Chapter 2 first discusses a number of basic concepts 
and results from probability theory that we will use. The most simple interesting queueing 
model is treated in chapter 4, and its multi server version is treated in the next chapter. 
Models with more general service or interarrival time distributions are analysed in the 
chapters 6, 7 and 8. Some simple variations on these models are discussed in chapter 10. 
Chapter 9 is devoted to queueing models with priority rules. The last chapter discusses 
some insentive systems. 

The text contains a lot of exercises and the reader is urged to try these exercises. This 
is really necessary to acquire skills to model and analyse new situations. 

1.1 Examples 

Below we briefly describe some situations in which queueing is important. 

Example 1.1.1 Supermarket. 

How long do customers have to wait at the checkouts? What happens with the waiting 



time during peak-hours? Are there enough checkouts? 

Example 1.1.2 Production system. 

A machine produces different types of products. 

What is the production lead time of an order? What is the reduction in the lead time 

when we have an extra machine? Should we assign priorities to the orders? 

Example 1.1.3 Post office. 

In a post office there are counters specialized in e.g. stamps, packages, financial transac- 
tions, etc. 

Are there enough counters? Separate queues or one common queue in front of counters 
with the same specialization? 

Example 1.1.4 Data communication. 

In computer communication networks standard packages called cells are transmitted over 

links from one switch to the next. In each switch incoming cells can be buffered when the 

incoming demand exceeds the link capacity. Once the buffer is full incoming cells will be 

lost. 

What is the cell delay at the switches? What is the fraction of cells that will be lost? What 

is a good size of the buffer? 

Example 1.1.5 Parking place. 

They are going to make a new parking place in front of a super market. 

How large should it be? 

Example 1.1.6 Assembly of printed circuit boards. 

Mounting vertical components on printed circuit boards is done in an assembly center 

consisting of a number of parallel insertion machines. Each machine has a magazine to 

store components. 

What is the production lead time of the printed circuit boards? How should the components 

necessary for the assembly of printed circuit boards be divided among the machines? 

Example 1.1.7 Call centers of an insurance company. 

Questions by phone, regarding insurance conditions, are handled by a call center. This call 
center has a team structure, where each team helps customers from a specific region only. 
How long do customers have to wait before an operator becomes available? Is the number 
of incoming telephone lines enough? Are there enough operators? Pooling teams? 

Example 1.1.8 Main frame computer. 

Many cashomats are connected to a big main frame computer handling all financial trans- 
actions. 

Is the capacity of the main frame computer sufficient? What happens when the use of 
cashomats increases? 
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Example 1.1.9 Toll booths. 

Motorists have to pay toll in order to pass a bridge. Are there enough toll booths? 

Example 1.1.10 Traffic lights. 

How do we have to regulate traffic lights such that the waiting times are acceptable? 
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Chapter 2 

Basic concepts from probability 
theory 

This chapter is devoted to some basic concepts from probability theory. 

2.1 Random variable 

Random variables are denoted by capitals, X , Y, etc. The expected value or mean of X is 
denoted by E(X) and its variance by o~ 2 (X) where cr(X) is the standard deviation of X. 
An important quantity is the coefficient of variation of the positive random variable X 
defined as 

a{X) 
ex 



E(xy 



The coefficient of variation is a (dimensionless) measure of the variability of the random 
variable X . 



2.2 Generating function 

Let X be a nonnegative discrete random variable with P(X = n) = p(n), n = 0, 1, 2, 
Then the generating function P x (z) of X is defined as 

oo 

P x {z) = E{z x ) = Y J P(n)z n . 

n=0 

Note that \P x (z)\ < 1 for all \z\ < 1. Further 

P x (0)=p(0), P x (l) = l, P X (1) = E(X), 
and, more general, 

P { x k \l) = E(X(X - 1) • • • (X - k + 1)), 
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where the superscript (k) denotes the kth derivative. For the generating function of the 
sum Z = X + Y of two independent discrete random variables X and Y, it holds that 

Pz(z) = Px(z)-Py(z). 
When Z is with probability q equal to X and with probability 1 — q equal to Y, then 

Pz(z) = qP X (z) + (l-q)P Y (z). 

2.3 Laplace- St ieltjes transform 

The Laplace-Stieltjes transform X(s) of a nonnegative random variable X with distribution 
function F(-), is defined as 

X(s) = E(e~ sX ) = / e- sx dF(x), s > 0. 

Jx=0 

When the random variable X has a density /(•), then the transform simplifies to 

fOO 

X(s) = / e~ sx f(x)dx, s > 0. 

Note that \X(s)\ < 1 for all s > 0. Further 

X(0) = 1, X'(0) = -E(X), X (k \0) = (-l) k E(X k ). 

For the transform of the sum Z = X + Y of two independent random variables X and Y, 
it holds that 

Z{s) = X{s)-Y{s). 

When Z is with probability q equal to X and with probability 1 — q equal to Y, then 

Z(s) = qX(s) + (l-q)Y(s). 

2.4 Useful probability distributions 

This section discusses a number of important distributions which have been found useful 
for describing random variables in many applications. 

2.4.1 Geometric distribution 

A geometric random variable X with parameter p has probability distribution 

p(X = n) = (1 -p)p n , n = 0,1,2,... 
For this distribution we have 

Px(z) = ^, E(X) = -*-, - 2 W = 7T J ^' c x = ~- 
1 — pz 1 — p (1 — p) 1 p 
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2.4.2 Poisson distribution 

A Poisson random variable X with parameter jj has probability distribution 

P(X = n) = t-re-t 1 , n = 0,l,2,... 

For the Poisson distribution it holds that 

P x (z) = e-rt-K E(X) = a\X) = p, c\ = -. 

2.4.3 Exponential distribution 

The density of an exponential distribution with parameter [x is given by 

f(t) = ixz-»\ t > 0. 
The distribution function equals 

F(t) = 1 - e _Ati , t > 0. 
For this distribution we have 

X(8) = -£-, E(X) = -, a 2 (X) = - 2 , c x = l. 

H + s fi fi z 

An important property of an exponential random variable X with parameter \x is the 
memoryless property. This property states that for all x > and t > 0, 



P(X > x + t\X > t) = P{X > x) = e 



■fix 



So the remaining lifetime of X, given that X is still alive at time t, is again exponentially 
distributed with the same mean 1/fi. We often use the memoryless property in the form 

P(X <t + At\X >t) = l-e- /iA ' = / uAt + o(At), (At->0), (2.1) 

where o(At), (At — > 0), is a shorthand notation for a function, g(At) say, for which 
g(At)/At tends to when At -> (see e.g. [ ]). 

If Xi, . . . , A n are independent exponential random variables with parameters ;Ui, . . . , yU n 
respectively, then min(Xi, . . . , X n ) is again an exponential random variable with parameter 
Ati + • • • + Un an d the probability that Xi is the smallest one is given by ^ij[}i\ -\ — • + /U n ), 
i = 1, . . . ,n. (see exercise 1). 
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2.4.4 Erlang distribution 

A random variable X has an Erlang-k (k = 1,2,...) distribution with mean k/fi if X 
is the sum of k independent random variables X±, . . . , Xj. having a common exponential 
distribution with mean 1/ pi. The common notation is E^n) or briefly E^. The density of 
an Ek(n) distribution is given by 

fa) = ^^9^, t > o. 

The distribution function equals 

F { t ) = l-Y J ^f e -^\ t>0. 

The parameter /x is called the scale parameter, k is the shape parameter. A phase diagram 
of the Ek distribution is shown in figure 2.1. 

-^ 1 1 J 2 1 ^ • • • U k 






Figure 2.1: Phase diagram for the Erlang- k distribution with scale parameter /i 

In figure 2.2 we display the density of the Erlang-fc distribution with mean 1 (so /x = k) 
for various values of k. 

The mean, variance and squared coefficient of variation are equal to 

e(x) = l, Ax) = * , 4 - \- 

The Laplace-Stieltjes transform is given by 

^ -(-£)■ 

A convenient distribution arises when we mix an E^_i and E^ distribution with the 
same scale parameters. The notation used is E^_\^. A random variable X has an E\._\j.{p) 
distribution, if X is with probability p (resp. 1 — p) the sum of k — 1 (resp. k) independent 
exponentials with common mean l//x The density of this distribution has the form 

where < p < 1. As p runs from 1 to 0, the squared coefficient of variation of the 
mixed Erlang distribution varies from l/(k — 1) to 1/k. It will appear (later on) that this 
distribution is useful for fitting a distribution if only the first two moments of a random 
variable are known. 
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0.5 1 1.5 2 2.5 3 3.5 4 

Figure 2.2: The density of the Erlang-A; distribution with mean 1 for various values of k 



2.4.5 Hyperexponential distribution 

A random variable X is hyperexponentially distributed if X is with probability pi, i = 
1, . . . , k an exponential random variable Xi with mean 1/yUj. For this random variable we 
use the notation Hk(pi, . . . ,Pk', jUi, • • • , ^k), or simply Hk- The density is given by 



/(*) = J2PiVie~ 



■flit 



t> 0, 



i=l 



and the mean is equal to 



i=l ri 

The Laplace-Stieltjes transform satisfies 



*(«) = £ 



PiAii 



i A*i + s 

The coefficient of variation ex of this distribution is always greater than or equal to 1 
(see exercise 3). A phase diagram of the Hk distribution is shown in figure 2.3. 
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Figure 2.3: Phase diagram for the hyperexponential distribution 

2.4.6 Phase-type distribution 

The preceding distributions are all special cases of the phase-type distribution. The notation 
is PH. This distribution is characterized by a Markov chain with states 1, . . . , k (the so- 
called phases) and a transition probability matrix P which is transient. This means that 
P n tends to zero as n tends to infinity. In words, eventually you will always leave the 
Markov chain. The residence time in state i is exponentially distributed with mean 1/yUj, 
and the Markov chain is entered with probability pi in state i, i = 1, . . . , k. Then the 
random variable X has a phase-type distribution if X is the total residence time in the 
preceding Markov chain, i.e. X is the total time elapsing from start in the Markov chain 
till departure from the Markov chain. 

We mention two important classes of phase-type distributions which are dense in the 
class of all non-negative distribution functions. This is meant in the sense that for any 
non-negative distribution function F(-) a sequence of phase- type distributions can be found 
which pointwise converges at the points of continuity of F(-). The denseness of the two 
classes makes them very useful as a practical modelling tool. A proof of the denseness can 
be found in [23, 24]. The first class is the class of Coxian distributions, notation Ck, and 
the other class consists of mixtures of Erlang distributions with the same scale parameters. 
The phase representations of these two classes are shown in the figures 2.4 and 2.5. 




1-Pi 



1-P2 



Pk-l 



l-Pk- 



\h 




Figure 2.4: Phase diagram for the Coxian distribution 

A random variable X has a Coxian distribution of order k if it has to go through up to 
at most k exponential phases. The mean length of phase n is 1//U n , n = 1, . . . , k. It starts 
in phase 1 . After phase n it comes to an end with probability \—p n and it enters the next 
phase with probability p n . Obviously pk = 0. For the Coxian-2 distribution it holds that 
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k 
n 

Figure 2.5: Phase diagram for the mixed Erlang distribution 

the squared coefficient of variation is greater than or equal to 0.5 (see exercise 8). 

A random variable X has a mixed Erlang distribution of order k if it is with probability 
p n the sum of n exponentials with the same mean 1//^, n = 1, . . . , k. 



2.5 Fitting distributions 

In practice it often occurs that the only information of random variables that is available 
is their mean and standard deviation, or if one is lucky, some real data. To obtain an 
approximating distribution it is common to fit a phase-type distribution on the mean, 
E(X), and the coefficient of variation, ex, of a given positive random variable X, by using 
the following simple approach. 

In case < ex < 1 one fits an Ek-x,k distribution (see subsection 2.4.4). More specifi- 
cally, if 

1 2 1 



k ~ x ~ Jfe-1' 

for certain k = 2,3, . . ., then the approximating distribution is with probability p (resp. 
1 — p) the sum of k — 1 (resp. k) independent exponentials with common mean l//i. By 
choosing (see e.g. [28]) 

P = y-^lkcx ~ W + 4) - k*c x Y% n = ^f y 

the Ek-i^k distribution matches E(X) and ex- 

In case Cx > 1 one fits a H 2 (pi,P2', ^1,^2) distribution. The hyperexponential distribu- 
tion however is not uniquely determined by its first two moments. In applications, the H2 
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distribution with balanced means is often used. This means that the normalization 
?1 = El 

is used. The parameters of the H 2 distribution with balanced means and fitting E(X) and 
c x (> 1) are given by 



pi= H 1t "n 



c 2 - 



2pi 2p 2 



£(X)' n E(X) 

In case c 2 ^ > 0.5 one can also use a Coxian-2 distribution for a two-moment fit. The 
following set is suggested by [18], 

Hi = 2/E(X), pi = 0.5/c|, // 2 = A*iPi- 

It also possible to make a more sophisticated use of phase-type distributions by, e.g., 
trying to match the first three (or even more) moments of X or to approximate the shape 
of X (see e.g. [29, 11, 13]). 

Phase-type distributions may of course also naturally arise in practical applications. 
For example, if the processing of a job involves performing several tasks, where each task 
takes an exponential amount of time, then the processing time can be described by an 
Erlang distribution. 

2.6 Poisson process 

Let N(t) be the number of arrivals in [0,t] for a Poisson process with rate A, i.e. the time 
between successive arrivals is exponentially distributed with parameter A and independent 
of the past. Then N(t) has a Poisson distribution with parameter At, so 

P(N(t) = k) = ( -^fe- xt , k = 0,1,2,... 

The mean, variance and coefficient of variation of N(t) are equal to (see subsection 2.4.2) 

1 



E(N(t)) = At, a 2 (N(t)) = At, c 2 N(t) - -. 

From (2.1) it is easily verified that 

P (arrival in (t, t + At]) = A At + o(At), (At -> 0). 

Hence, for small At, 

P(arrival in (t, t + At]) w AAt. (2.2) 
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So in each small time interval of length At the occurence of an arrival is equally likely. In 
other words, Poisson arrivals occur completely random in time. In figure 2.6 we show a 
realization of a Poisson process and an arrival process with Erlang-10 interarrival times. 
Both processes have rate 1. The figure illustrates that Erlang arrivals are much more 
equally spread out over time than Poisson arrivals. 

P i SS0n I I I 111} U I LU LU I U 



Erlang-10 





t 



Figure 2.6: A realization of Poisson arrivals and Erlang-10 arrivals, both with rate 1 

The Poisson process is an extremely useful process for modelling purposes in many 
practical applications, such as, e.g. to model arrival processes for queueing models or 
demand processes for inventory systems. It is empirically found that in many circumstances 
the arising stochastic processes can be well approximated by a Poisson process. 

Next we mention two important properties of a Poisson process (see e.g. [ ]). 

(i) Merging. 

Suppose that N\(t) and N 2 (t) are two independent Poisson processes with respective 
rates Ai and A2. Then the sum Ni(t) + AT 2 (t) of the two processes is again a Poisson 
process with rate Ai + A2. 

(ii) Splitting. 

Suppose that N(t) is a Poisson process with rate A and that each arrival is marked 
with probability p independent of all other arrivals. Let Ni(t) and A^t) denote 
respectively the number of marked and unmarked arrivals in [0, £]. Then N\(t) and 
N 2 (t) are both Poisson processes with respective rates Xp and A(l — p). And these 
two processes are independent. 

So Poisson processes remain Poisson processes under merging and splitting. 
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2.7 Exercises 

Exercise 1. 

Let Xi, . . . ,X n be independent exponential random variables with mean E(Xj) = 1//Uj, 
i = 1, . . . ,n. Define 

Y n = min(Xi, . . . , X n ), Z n = max(Xi, . . . , X n ). 

(i) Determine the distributions of Y n and Z n . 

(ii) Show that the probability that Xi is the smallest one among X\, . . . , X n is equal to 
fJ-i/ifJ-i -\ h/i n ), i = l,...,n. 

Exercise 2. 

Let X\, X2, ... be independent exponential random variables with mean 1/jj, and let iV be 

a discrete random variable with 

P(Ar = fc ) = (l-p)p fc - 1 , fc = l,2,..., 

where < p < 1 (i.e. A^ is a shifted geometric random variable). Show that S defined as 

TV 
n=l 

is again exponentially distributed with parameter (1 —p)/J,. 

Exercise 3. 

Show that the coefficient of variation of a hyperexponential distribution is greater than or 

equal to 1. 

Exercise 4. (Poisson process) 

Suppose that arrivals occur at 21, T2, . . .. The interarrival times A n = T n — T n _i are 
independent and have common exponential distribution with mean 1/A, where To = by 
convention. Let N(t) denote the number of arrivals is [0, t] and define for n = 0, 1, 2, . . . 

p n (t) = P(N(t) = n), t> 0. 

(i) Determine Po{t)- 
(ii) Show that for n = 1, 2, . . . 

P ' n (t) = -X Pn (t) + Xpn.^t), t>0, 

with initial condition p n (0) = 0. 
(iii) Solve the preceding differential equations for n = 1, 2, . . . 
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Exercise 5. (Poisson process) 

Suppose that arrivals occur at 71, T2, . . .. The interarrival times A n = T n — T n _i are 
independent and have common exponential distribution with mean 1/A, where To = by 
convention. Let N(t) denote the number of arrivals is [0, t] and define for n = 0, 1, 2, . . . 

Pn (t) = P(N(t) = n), t> 0. 

(i) Determine Po(t). 
(ii) Show that for n = 1, 2, . . . 

rt 

Pn(t) = / p n -i(t - x)\e~ Xx dx, t > 0. 

JO 

(iii) Solve the preceding integral equations for n = 1, 2, . . . 

Exercise 6. (Poisson process) 

Prove the properties (i) and (ii) of Poisson processes, formulated in section 2.6. 

Exercise 7. (Fitting a distribution) 

Suppose that processing a job on a certain machine takes on the average 4 minutes with a 
standard deviation of 3 minutes. Show that if we model the processing time as a mixture 
of an Erlang-1 (exponential) distribution and an Erlang-2 distribution with density 

/(i)=p/ie-^ + (l-p)/i 2 te-^ 

the parameters p and \i can be chosen in such a way that this distribution matches the 
mean and standard deviation of the processing times on the machine. 

Exercise 8. 

Consider a random variable X with a Coxian-2 distribution with parameters /xi and //2 

and branching probability p\ . 

(i) Show that c\ > 0.5. 

(ii) Show that if /Ji < fi2, then this Coxian-2 distribution is identical to the Coxian- 
2 distribution with parameters fa, fa and f>\ where fa = // 2 , fa = fJ-i and p\ = 
1 - (l-pi)/!i//i 2 . 

Part (ii) implies that for any Coxian-2 distribution we may assume without loss of generality 

that yUi > yU2- 

Exercise 9. 

Let X and Y be exponentials with parameters \x and A, respectively. Suppose that A < [x. 
Let Z be equal to X with probability X//j and equal to X + Y with probability 1 — A/ fi. 
Show that Z is an exponential with parameter A. 
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Exercise 10. 

Consider a Hi distribution with parameters \X\ > \x% and branching probabilities q\ and 
q2, respectively. Show that the C 2 distribution with parameters \x\ and \x^ and branching 
probability px given by 

p x = 1 - (gi/^i + q2ixi)l\xi , 

is equivalent to the H 2 distribution. 

Exercise 11. (Poisson distribution) 

Let X\, . . . , X n be independent Poisson random variables with means //1, . . . , // n , respec- 
tively. Show that the sum X\ + ■ ■ ■ + X n is Poisson distributed with mean n^ + • • • + /i ra . 
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Chapter 3 

Queueing models and some 
fundamental relations 



In this chapter we describe the basic queueing model and we discuss some important fun- 
damental relations for this model. These results can be found in every standard textbook 
on this topic, see e.g. [14, 20, 28]. 



3.1 Queueing models and Kendall's notation 

The basic queueing model is shown in figure 3.1. It can be used to model, e.g., machines 
or operators processing orders or communication equipment processing information. 




Figure 3.1: Basic queueing model 



Among others, a queueing model is characterized by: 

• The arrival process of customers. 

Usually we assume that the interarrival times are independent and have a common 
distribution. In many practical situations customers arrive according to a Poisson 
stream (i.e. exponential interarrival times). Customers may arrive one by one, or 
in batches. An example of batch arrivals is the customs office at the border where 
travel documents of bus passengers have to be checked. 
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• The behaviour of customers. 

Customers may be patient and willing to wait (for a long time). Or customers may 
be impatient and leave after a while. For example, in call centers, customers will 
hang up when they have to wait too long before an operator is available, and they 
possibly try again after a while. 

• The service times. 

Usually we assume that the service times are independent and identically distributed, 
and that they are independent of the interarrival times. For example, the service 
times can be deterministic or exponentially distributed. It can also occur that service 
times are dependent of the queue length. For example, the processing rates of the 
machines in a production system can be increased once the number of jobs waiting 
to be processed becomes too large. 

• The service discipline. 

Customers can be served one by one or in batches. We have many possibilities for 
the order in which they enter service. We mention: 

— first come first served, i.e. in order of arrival; 

— random order; 

— last come first served (e.g. in a computer stack or a shunt buffer in a production 
line); 

— priorities (e.g. rush orders first, shortest processing time first); 

— processor sharing (in computers that equally divide their processing power over 
all jobs in the system). 

• The service capacity. 

There may be a single server or a group of servers helping the customers. 

• The waiting room. 

There can be limitations with respect to the number of customers in the system. For 
example, in a data communication network, only finitely many cells can be buffered 
in a switch. The determination of good buffer sizes is an important issue in the design 
of these networks. 

Kendall introduced a shorthand notation to characterize a range of these queueing mod- 
els. It is a three-part code a/b/c. The first letter specifies the interarrival time distribution 
and the second one the service time distribution. For example, for a general distribution 
the letter G is used, M for the exponential distribution (M stands for Memoryless) and 
D for deterministic times. The third and last letter specifies the number of servers. Some 
examples are M/M/l, M/M/c, M/G/l, G/M/l and M/D/l. The notation can be ex- 
tended with an extra letter to cover other queueing models. For example, a system with 
exponential interarrival and service times, one server and having waiting room only for N 
customers (including the one in service) is abbreviated by the four letter code M/M/l/N. 
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In the basic model, customers arrive one by one and they are always allowed to enter 
the system, there is always room, there are no priority rules and customers are served in 
order of arrival. It will be explicitly indicated (e.g. by additional letters) when one of these 
assumptions does not hold. 

3.2 Occupation rate 

In a single-server system G/G/l with arrival rate A and mean service time E(B) the 
amount of work arriving per unit time equals XE(B). The server can handle 1 unit work 
per unit time. To avoid that the queue eventually grows to infinity, we have to require that 
XE(B) < 1. Without going into details, we note that the mean queue length also explodes 
when XE(B) = 1, except in the D/D/l system, i.e., the system with no randomness at all. 
It is common to use the notation 

p = XE(B). 

If p < 1, then p is called the occupation rate or server utilization, because it is the fraction 
of time the server is working. 

In a multi-server system G/G/c we have to require that XE(B) < c. Here the occupa- 
tion rate per server is p = XE(B)/c. 

3.3 Performance measures 

Relevant performance measures in the analysis of queueing models are: 

• The distribution of the waiting time and the sojourn time of a customer. The sojourn 
time is the waiting time plus the service time. 

• The distribution of the number of customers in the system (including or excluding 
the one or those in service). 

• The distribution of the amount of work in the system. That is the sum of service times 
of the waiting customers and the residual service time of the customer in service. 

• The distribution of the busy period of the server. This is a period of time during 
which the server is working continuously. 

In particular, we are interested in mean performance measures, such as the mean waiting 
time and the mean sojourn time. 

Now consider the G/G/c queue. Let the random variable L(t) denote the number of 
customers in the system at time i, and let S n denote the sojourn time of the nth customer 
in the system. Under the assumption that the occupation rate per server is less than one, 
it can be shown that these random variables have a limiting distribution as t — > oo and 
n — ► oo. These distributions are independent of the initial condition of the system. 
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Let the random variables L and S have the limiting distributions of L(t) and 5 n , 
respectively. So 

p k = P(L = k) = lim P(L(t) = k), F s (x) = P(S <x)= lim P(S n < x). 

t—>oo n— >oo 

The probability pk can be interpreted as the fraction of time that k customers are in the 
system, and Fs(x) gives the probability that the sojourn time of an arbitrary customer 
entering the system is not greater than x units of time. It further holds with probability 
1 that 



I r t i n 

lim - / L(x)dx = E(L), lim - V S k = E(S). 

t-oo t Jx=o v ' y ' ^njg 



So the long-run average number of customers in the system and the long-run average 
sojourn time are equal to E(L) and E(S), respectively. A very useful result for queueing 
systems relating E(L) and E(S) is presented in the following section. 

3.4 Little's law 

Little's law gives a very important relation between E(L), the mean number of customers 
in the system, E(S), the mean sojourn time and A, the average number of customers 
entering the system per unit time. Little's law states that 

E(L) = XE(S). (3.1) 

Here it is assumed that the capacity of the system is sufficient to deal with the customers 
(i.e. the number of customers in the system does not grow to infinity). 

Intuitively, this result can be understood as follows. Suppose that all customers pay 1 
dollar per unit time while in the system. This money can be earned in two ways. The first 
possibility is to let pay all customers "continuously" in time. Then the average reward 
earned by the system equals E(L) dollar per unit time. The second possibility is to let 
customers pay 1 dollar per unit time for their residence in the system when they leave. In 
equilibrium, the average number of customers leaving the system per unit time is equal 
to the average number of customers entering the system. So the system earns an average 
reward of XE(S) dollar per unit time. Obviously, the system earns the same in both cases. 
For a rigorous proof, see [17, 25]. 

To demonstrate the use of Little's law we consider the basic queueing model in figure 
3.1 with one server. For this model we can derive relations between several performance 
measures by applying Little's law to suitably defined (sub)systems. Application of Little's 
law to the system consisting of queue plus server yields relation (3.1). Applying Little's 
law to the queue (excluding the server) yields a relation between the queue length L q and 
the waiting time W, namely 

E{L q ) = XE(W). 
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Finally, when we apply Little's law to the server only, we obtain (cf. section 3.2) 

p = \E(B), 

where p is the mean number of customers at the server (which is the same as the fraction 
of time the server is working) and E(B) the mean service time. 

3.5 PASTA property 

For queueing systems with Poisson arrivals, so for M / ■ / ■ systems, the very special property 
holds that arriving customers find on average the same situation in the queueing system 
as an outside observer looking at the system at an arbitrary point in time. More precisely, 
the fraction of customers finding on arrival the system in some state A is exactly the same 
as the fraction of time the system is in state A. This property is only true for Poisson 
arrivals. 

In general this property is not true. For instance, in a D/D/l system which is empty 
at time 0, and with arrivals at 1, 3, 5, . . . and service times 1, every arriving customer finds 
an empty system, whereas the fraction of time the system is empty is 1/2. 

This property of Poisson arrivals is called PASTA property, which is the acrynom for 
Poisson Arrivals See Time Averages. Intuitively, this property can be explained by the 
fact that Poisson arrivals occur completely random in time (see (2.2)). A rigorous proof of 
the PASTA property can be found in [31, 32]. 

In the following chapters we will show that in many queueing models it is possible to 
determine mean performance measures, such as E(S) and E(L), directly (i.e. not from 
the distribution of these measures) by using the PASTA property and Little's law. This 
powerful approach is called the mean value approach. 
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3.6 Exercises 

Exercise 12. 

In a gas station there is one gas pump. Cars arrive at the gas station according to a Poisson 
proces. The arrival rate is 20 cars per hour. An arriving car finding n cars at the station 
immediately leaves with probability q n = n/4, and joins the queue with probability 1 — q n , 
n = 0, 1, 2, 3,4. Cars are served in order of arrival. The service time (i.e. the time needed 
for pumping and paying) is exponential. The mean service time is 3 minutes. 

(i) Determine the stationary distribution of the number of cars at the gas station. 

(ii) Determine the mean number of cars at the gas station. 

(iii) Determine the mean sojourn time (waiting time plus service time) of cars deciding 
to take gas at the station. 

(iv) Determine the mean sojourn time and the mean waiting time of all cars arriving at 
the gas station. 
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Chapter 4 
M/M/l queue 



In this chapter we will analyze the model with exponential interarrival times with mean 
1/A, exponential service times with mean 1/p and a single server. Customers are served 
in order of arrival. We require that 

A r 
p=-<l, 

since, otherwise, the queue length will explode (see section 3.2). The quantity p is the 
fraction of time the server is working. In the following section we will first study the 
time- dependent behaviour of this system. After that, we consider the limiting behaviour. 

4.1 Time-dependent behaviour 

The exponential distribution allows for a very simple description of the state of the system 
at time t, namely the number of customers in the system (i.e. the customers waiting in 
the queue and the one being served). Neither we do have to remember when the last 
customer arrived nor we have to register when the last customer entered service. Since the 
exponential distribution is memoryless (see 2.1), this information does not yield a better 
prediction of the future. 

Let p n (t) denote the probability that at time t there are n customers in the system, 
n = 0, 1, . . . Based on property (2.1) we get, for At — > 0, 

p (t + At) = (l-XAt) P o(t)+pAt Pl (t) + o(At), 
Pn {t + At) = XAtp n _ l (t) + (l-{X + p)At)p n {t)+pAt Pn+l (t) + o{At), 
n = l,2,... 

Hence, by letting At — ► 0, we obtain the following infinite set of differential equations for 
the probabilities p n (t). 

Po(t) = -Apo(*) + m(*), / 41 \ 

P' n (t) = Ap n _ 1 (t)-(A + /i) Pn (t)+/ip n+1 (t), n=l,2,... { ' ] 
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It is difficult to solve these differential equations. An explicit solution for the probabilities 
p n {t) can be found in [11] (see p. 77). The expression presented there is an infinite sum 
of modified Bessel functions. So already one of the simplest interesting queueing models 
leads to a difficult expression for the time-dependent behavior of its state probabilities. For 
more general systems we can only expect more complexity. Therefore, in the remainder 
we will focus on the limiting or equilibrium behavior of this system, which appears to be 
much easier to analyse. 



4.2 Limiting behavior 

One may show that as t —> oo, then p' n (t) — > and p n (t) — > p n (see e.g. [ ■]). Hence, from 
(4.1) it follows that the limiting or equilibrium probabilities p n satisfy the equations 

= -Xpo + /ipi, (4.2) 

= Xpn-l - (A + fj)p n + (J>Pn+U n = 1, 2, . . . (4.3) 

Clearly, the probabilities p n also satisfy 



I>n=l, (4-4) 

71=0 

which is called the normalization equation. It is also possible to derive the equations (4.2) 
and (4.3) directly from a flow diagram, as shown in figure 4.1. 

,wo ... (5Q3 ... 

a a a 

Figure 4.1: Flow diagram for the M/M/l model 

The arrows indicate possible transitions. The rate at which a transition occurs is A for 
a transition from nton + 1 (an arrival) and [x for a transition from n + 1 to n (a departure). 
The number of transitions per unit time from nton + 1, which is also called the flow from 
n to n + 1, is equal to p n , the fraction of time the system is in state n, times A, the rate 
at arrivals occur while the system is in state n. The equilibrium equations (4.2) and (4.3) 
follow by equating the flow out of state n and the flow into state n. 

For this simple model there are many ways to determine the solution of the equations 
(4.2)-(4.4). Below we discuss several approaches. 
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4.2.1 Direct approach 

The equations (4.3) are a second order recurrence relation with constant coefficients. Its 
general solution is of the form 

p n = cix" + c 2 X2 , n = 0,1,2,... (4.5) 

where x\ and x 2 are roots of the quadratic equation 

A - (A + p)x + fix 2 = 0. 

This equation has two zeros, namely x = 1 and x = X/fx = p. So all solutions to (4.3) are 
of the form 

Pn = C 1 +C 2 p n , 72 = 0,1,2, ... 

Equation (4.4), stating that the sum of all probabilities is equal to 1, of course directly 
implies that C\ must be equal to 0. That C\ must be equal to also follows from (4.2) by 
substituting the solution (4.5) into (4.2). 

The coefficient c 2 finally follows from the normalization equation (4.4), yielding that 
c 2 = 1 — p. So we can conclude that 

Pn = (l-p)p n , n = 0,1,2,... (4.6) 

Apparantly, the equilibrium distribution depends upon A and p only through their ratio p. 

4.2.2 Recursion 

One can use (4.2) to express p\ in p yielding 

Pi = PPo- 
Substitution of this relation into (4.3) for n = 1 gives 

V2 = p 2 Po- 

By substituting the relations above into (4.3) for n = 2 we obtain ps, and so on. Hence we 
can recursively express all probabilities in terms of p , yielding 

Pn = p n Po, n = 0,1,2, .... 

The probability po finally follows from the normalization equation (4.4). 
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4.2.3 Generating function approach 

The probability generating function of the random variable L, the number of customers in 
the system, is given by 

oo 

P L {z) = Y,PnZ n , (4.7) 

ra=0 

which is properly defined for z with \z\ < 1. By multiplying the nth equilibrium equation 
with z n and then summing the equations over all n, the equilibrium equations for p n can 
be transformed into the following single equation for Pl{z), 

= fipo(l - z- 1 ) + (Xz + pz~ x - (A + n))P L {z). 

The solution of this equation is 

Pl(z) = r ^- = I—?- = £ (1 - p)p n z\ (4.8) 

1-pz 1-pz n=0 

where we used that P(l) = 1 to determine p = 1 — p (cf. section 3.2). Hence, by equating 
the coefficients of z n in (4.7) and (4.8) we retrieve the solution (4.6). 

4.2.4 Global balance principle 

The global balance principle states that for each set of states A, the flow out of set A is 
equal to the flow into that set. In fact, the equilibrium equations (4.2)-(4.3) follow by 
applying this principle to a single state. But if we apply the balance principle to the set 
A = {0,l,...,n — 1} we get the very simple relation 

Xp n -l = pp n , H = 1,2,... 

Repeated application of this relation yields 

Pn = p n po, n = 0,1,2,... 

so that, after normalization, the solution (4.6) follows. 

4.3 Mean performance measures 

From the equilibrium probabilities we can derive expressions for the mean number of cus- 
tomers in the system and the mean time spent in the system. For the first one we get 



P_ 
>-9 



E ( L ) = J2 n Pn = Y 



and by applying Little's law, 

E(S) = -^-. (4.9) 

1 - p 
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If we look at the expressions for E(L) and E(S) we see that both quantities grow to infinity 
as p approaches unity. The dramatic behavior is caused by the variation in the arrival and 
service process. This type of behavior with respect to p is characteristic for almost every 
queueing system. 

In fact, E(L) and E(S) can also be determined directly, i.e. without knowing the 
probabilities p n , by combining Little's law and the PASTA property (see section 3.5). 
Based on PASTA we know that the average number of customers in the system seen by an 
arriving customer equals E(L) and each of them (also the one in service) has a (residual) 
service time with mean 1/p. The customer further has to wait for its own service time. 
Hence 

E(S) = E(L)- + - . 
p p 

This relation is known as the arrival relation. Together with 

E(L) = XE(S) 

we find expression (4.9). This approach is called the mean value approach. 

The mean number of customers in the queue, E(L q ), can be obtained from E(L) by 
subtracting the mean number of customers in service, so 

E{L*) = E{L)-p=-?—. 

\- p 

The mean waiting time, E(W), follows from E(S) by subtracting the mean service time 
(or from E(L q ) by applying Little's law). This yields 

pin 



E(W) = E{S) - 1/p 



1-p 



4.4 Distribution of the sojourn time and the waiting 
time 

It is also possible to derive the distribution of the sojourn time. Denote by L a the number of 
customers in the system just before the arrival of a customer and let Bk be the service time 
of the fcth customer. Of course, the customer in service has a residual service time instead 
of an ordinary service time. But these are the same, since the exponential service time 
distribution is memoryless. So the random variables B k are independent and exponentially 
distributed with mean 1/p. Then we have 

S = £ B k . (4.10) 

fe=i 

By conditioning on L a and using that L a and Bk are independent it follows that 

L a +1 oo n+1 

p(s > t) = p( E B k >t) = Y, p(E B * > *) p ( L ° = n )- ( 4 - n ) 

fc=l n=0 k=l 
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The problem is to find the probability that an arriving customer finds n customers in the 
system. PASTA states that the fraction of customers finding on arrival n customers in the 
system is equal to the fraction of time there are n customers in the system, so 

P(L a = n)=p n = (l-p)p n . (4.12) 

Substituting (4.12) in (4.11) and using that YJl=\B k is Erlang-(n + 1) distributed, yields 
(cf. exercise 2) 

P(S>t) = EE^e-^i-pK 

n=0 fc=0 K - 

oo oo I ,\k 

= EE^(i-p) P " 

k=0n=k "" 



E 






k=0 

= e - Kl - p) \ t > 0. (4.13) 

Hence, S is exponentially distributed with parameter /i(l — p). This result can also be 
obtained via the use of transforms. From (4.10) it follows, by conditioning on L a , that 

S(s) = E(e~ sS ) 

oo 

= E P ( La = n)E{e- s(Bl+ - +Bn+l) ) 

n=0 

oo 

= Y,( 1 -p)p nE ( e ~ sBl )--- E ( e ~ sBn+1 )- 

n=0 

Since B k is exponentially distributed with parameter fi, we have (see subsection 2.4.3) 

sB k , _ P 



SO 



E(e- SD ») - 

P 



oo / \ n +l /I \ 

s(s) = y.^-p)p J " ^ " ( P) 



n=0 



JX + S J /i(l - p) + s ' 



from which we can conclude that 5" is an exponential random variable with parameter 
/i(l-p). 

To find the distribution of the waiting time W, note that S = W + B, where the random 
variable B is the service time. Since W and B are independent, it follows that 



S(s) = 


--W(s) 


■ B(s) - 


= W{s) ■ 


A* 

P + 


S 






and thus, 
















W(s) ■■ 


(1- 


P)(P + 


!U 


~P)- 


1 + P- 


p(l 
p(l- 


~P) 


n(i 


~P) + 


p) + s 
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From the transform of W we conclude (see subsection 2.3) that W is with probability 
(1 — p) equal to zero, and with probability p equal to an exponential random variable with 
parameter p(l — p). Hence 

P{W >t) = pe'^ 1 -^, t > 0. (4.14) 

The distribution of W can, of course, also be obtained along the same lines as (4.13). 
Note that 

p(w > t\w > o) = gjj^l = „-*-*, 

so the conditional waiting time W\W > is exponentially distributed with parameter 
/i(l-p). 

In table 4.1 we list for increasing values of p the mean waiting time and some waiting 
time probabilities. From these results we see that randomness in the arrival and service 
process leads to (long) waiting times and the waiting times explode as the server utilization 
tends to one. 



p 


E{W) 


P(W > t) 




t 5 10 20 


0.5 


1 


0.04 0.00 0.00 


0.8 


4 


0.29 0.11 0.02 


0.9 


9 


0.55 0.33 0.12 


0.95 


19 


0.74 0.58 0.35 



Table 4.1: Performance characteristics for the M/M/l with mean service time 1 



Remark 4.4.1 (PASTA property) 

For the present model we can also derive relation (4.12) directly from the flow diagram 
4.1. Namely, the average number of customers per unit time finding on arrival n customers 
in the system is equal to Xp n . Dividing this number by the average number of customers 
arriving per unit time gives the desired fraction, so 

P(L a = n) = ^= Pn . 

4.5 Priorities 

In this section we consider an M/M/l system serving different types of customers. To keep 
it simple we suppose that there are two types only, type 1 and 2 say, but the analysis can 
easily be extended the situation with more types of customers (see also chapter 9). Type 1 
and type 2 customers arrive according to independent Poisson processes with rate Ai, and 
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A2 respectively. The service times of all customers are exponentially distributed with the 
same mean 1///. We assume that 

P1 + P2 < 1, 

where pi = \i/p, i.e. the occupation rate due to type i customers. Type 1 customers are 
treated with priority over type 2 jobs. In the following subsections we will consider two 
priority rules, preemptive-resume priority and non-preemptive priority. 

4.5.1 Preemptive-resume priority 

In the preemptive resume priority rule, type 1 customers have absolute priority over type 
2 jobs. Absolute priority means that when a type 2 customer is in service and a type 1 
customer arrives, the type 2 service is interrupted and the server proceeds with the type 1 
customer. Once there are no more type 1 customers in the system, the server resumes the 
service of the type 2 customer at the point where it was interrupted. 

Let the random variable Li denote the number of type i customers in the system and 
Si the sojourn time of a type i customer. Below we will determine E(Li) and E(Si) for 
% = 1,2. 

For type 1 customers the type 2 customers do not exist. Hence we immediately have 

E(S 1 ) = -^-, E(L 1 ) = -^-. (4.15) 

1 - Pi 1 - Pi 

Since the (residual) service times of all customers are exponentially distributed with the 
same mean, the total number of customers in the system does not depend on the order in 
which the customers are served. So this number is the same as in the system where all 
customers are served in order of arrival. Hence, 

E{L X ) + E(L 2 ) = - Pl+P2 , (4.16) 

1 - Pi - Pi 

and thus, inserting (4.15), 

Pi + P2 Pi P2 



E(L 2 



I-P1-P2 1-Pi (1 -pi)(l -pi -p 2 ) ' 
and applying Little's law, 

F(q v = MM = i/p 

[ 2) A 2 (i- pi )(i- pi - p2 )- 

Example 4.5.1 For Ai = 0.2, A 2 = 0.6 and \x = 1, we find in case all customers are 
treated in order of arrival, 

E(S) = — - — = 5, 
v ; 1-0.8 

and in case type 1 customers have absolute priority over type 2 jobs, 

E(Si) = = 1.25, E(S 3 ) = -, T7 r = 6.25. 

v ; 1-0.2 ' y ' (l-0.2)(l-0.8) 
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4.5.2 Non- preemptive priority 

We now consider the situation that type 1 customers have nearly absolute priority over 
type 2 jobs. The difference with the previous rule is that type 1 customers are not allowed 
to interrupt the service of a type 2 customers. This priority rule is therefore called non- 
preemptive. 

For the mean sojourn time of type 1 customers we find 

E(Si) = E(L 1 )- + - + p 2 -. 
(X ix ix 

The last term reflects that when an arriving type 1 customer finds a type 2 customer 
in service, he has to wait until the service of this type 2 customer has been completed. 
According to PASTA the probability that he finds a type 2 customer in service is equal 
to the fraction of time the server spends on type 2 customers, which is p 2 . Together with 
Little's law, 

£(Li) = \ X E(S X ), 
we obtain 

E( Sl )= il + p2)/ ^ E( Ll )= {1 + p2)pl . 
1 - Pi 1 - Pi 

For type 2 customers it follows from (4.16) that 

rpfr N (1-Pl(l-Pl-P 2 ))p2 

(1 -pi)(l -pi -p 2 ) 
and applying Little's law, 

evc n _ (1-Pi(l-Pi-P 2 ))/P 
E{S2) ~ (1-Pi)(l-Pi-P 2 ) • 

Example 4.5.2 For Ai = 0.2, A 2 = 0.6 and \x = 1, we get 

B(Sl )_I±™ 2, BW 1-0.2(1-0.8) . 



0.2 ' v ' (l-0.2)(l-0.8) 

4.6 Busy period 

In a servers life we can distinguish cycles. A cycle is the time that elapses between two 
consecutive arrivals finding an empty system. Clearly, a cycle starts with a busy period BP 
during which the server is helping customers, followed by an idle period IP during which 
the system is empty. 

Due to the memoryless property of the exponential distribution (see subsection 2.4.3), 
an idle period IP is exponentially distributed with mean 1/A. In the following subsections 
we determine the mean and the distribution of a busy period BP. 
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4.6.1 Mean busy period 

It is clear that the mean busy period divided by the mean cycle length is equal to the 
fraction of time the server is working, so 

E(BP) E(BP) 

P- 



E(BP) + E(IP) E(BP) + 1/X 
Hence, 



E(BP) 



P 



4.6.2 Distribution of the busy period 

Let the random variable C n be the time till the system is empty again if there are now n 
customers present in the system. Clearly, G\ is the length of a busy period, since a busy 
period starts when the first customer after an idle period arrives and it ends when the 
system is empty again. The random variables C n satisfy the following recursion relation. 
Suppose there are n(> 0) customers in the system. Then the next event occurs after an 
exponential time with parameter X + fj,: with probability X/(X + fx) a new customer arrives, 
and with probability ///(A + n) service is completed and a customer leaves the system. 
Hence, for n = 1, 2, . . ., 

q ^ x { Cn+i witn probability A/ (A + //), (A 17) 

n 1 C n -\ with probability p/(X + //), 

where X is an exponential random variable with parameter A + \x. From this relation we 
get for the Laplace-Stieltjes transform C n (s) of C n that 



A + /i + s \ A + /i X + n) 

and thus, after rewriting, 

(A + fi + s)C n (s) = XC n+1 (s) + /i(7 n _i(s), n = 1, 2, . . . 
For fixed s this equation is a very similar to (4.3). Its general solution is 

C n (s) = ciajj (s) + c 2 X2 (s), n = 0, 1, 2, . . . 
where Xi(s) and £2(5) are the roots of the quadratic equation 

(A + [i + s)x = Ax 2 + yu, 

satisfying < X\{s) < 1 < x-i{s)- Since < C n (s) < 1 it follows that c^ = 0. The 
coefficient C\ follows from the fact that Cq = and hence Cq(s) = 1, yielding c\ = 1. 
Hence we obtain 

£ n (a) = a;?(s), 
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and in particular, for the Laplace-Stieltjes transform BP(s) of the busy period BP, we 
find 

1 



BP(s) = Ci(s) = xi(s) = — (X + p + s-yJ(X + p + s) 2 - 4A/z) . 
By inverting this transform (see e.g. [■]) we get for the density fspit) of BP, 

f BP (t) = -^- e -( A+ ^ t / 1 (2t v / A^), t > 0, 
where I%(-) denotes the modified Bessel function of the first kind of order one, i.e. 

u (x) = T { ' ' 

In table 4.2 we list for some values of p the probability P(BP > t) for a number of t 
values. If you think of the situation that 1/// is one hour, then 10% of the busy periods 
lasts longer than 2 days (16 hours) and 5% percent even longer than 1 week, when p = 0.9. 
Since the mean busy period is 10 hours in this case, it is not unlikely that in a month time 
a busy period longer than a week occurs. 



p 






P(BP> 


*) 






1 


t 1 


2 


4 8 


16 


40 


80 


0.8 


0.50 


0.34 


0.22 0.13 


0.07 


0.02 


0.01 


0.9 


0.51 


0.36 


0.25 0.16 


0.10 


0.05 


0.03 


0.95 


0.52 


0.37 


0.26 0.18 


0.12 


0.07 


0.04 



Table 4.2: Probabilities for the busy period duration for the M/M/l with mean service 
time equal to 1 



4.7 Java applet 

For the performance evaluation of the M/M/l queue a JAVA applet is avalaible on the 
World Wide Web. The link to this applet is http://www.win.tue.nl/cow/Q2. The applet 
can be used to evaluate the mean value as well as the distribution of, e.g., the waiting time 
and the number of customers in the system. 
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4.8 Exercises 

Exercise 13. (bulk arrivals) 

In a work station orders arrive according to a Poisson arrival process with arrival rate A. 

An order consists of TV independent jobs. The distribution of TV is given by 



P(N = k) = {l-p)p k ~ 1 

with k = 1,2,... and < p < 1. Each job requires an exponentially distributed amount 
of processing time with mean 1/fi. 

(i) Derive the distribution of the total processing time of an order. 

(ii) Determine the distribution of the number of orders in the system. 

Exercise 14. (variable production rate) 

Consider a work station where jobs arrive according to a Poisson process with arrival rate 

A. The jobs have an exponentially distributed service time with mean 1/ ' \x. So the service 

completion rate (the rate at which jobs depart from the system) is equal to [x. 

If the queue length drops below the threshold Qi the service completion rate is lowered to 

Hl. If the queue length reaches Qh, where Qh > Ql, the service rate is increased to n H . 

(L stands for low, H for high.) 

Determine the queue length distribution and the mean time spent in the system. 

Exercise 15. 

A repair man fixes broken televisions. The repair time is exponentially distributed with 
a mean of 30 minutes. Broken televisions arrive at his repair shop according to a Poisson 
stream, on average 10 broken televisions per day (8 hours). 

(i) What is the fraction of time that the repair man has no work to do? 

(ii) How many televisions are, on average, at his repair shop? 

(iii) What is the mean throughput time (waiting time plus repair time) of a television? 

Exercise 16. 

In a gas station there is one gas pump. Cars arrive at the gas station according to a Poisson 
process. The arrival rate is 20 cars per hour. Cars are served in order of arrival. The service 
time (i.e. the time needed for pumping and paying) is exponentially distributed. The mean 
service time is 2 minutes. 

(i) Determine the distribution, mean and variance of the number of cars at the gas 
station. 

(ii) Determine the distribution of the sojourn time and the waiting time. 

(iii) What is the fraction of cars that has to wait longer than 2 minutes? 
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An arriving car finding 2 cars at the station immediately leaves. 

(iv) Determine the distribution, mean and variance of the number of cars at the gas 
station. 

(v) Determine the mean sojourn time and the mean waiting time of all cars (including 
the ones that immediately leave the gas station). 

Exercise 17. 

A gas station has two pumps, one for gas and the other for LPG. For each pump customers 
arrive according to a Poisson proces. On average 20 customers per hour for gas and 5 
customers for LPG. The service times are exponential. For both pumps the mean service 
time is 2 minutes. 

(i) Determine the distribution of the number of customers at the gas pump, and at the 
LPG pump. 

(ii) Determine the distribution of the total number of customers at the gas station. 

Exercise 18. 

Consider an M/M/l queue with two types of customers. The mean service time of all 
customers is 5 minutes. The arrival rate of type 1 customers is 4 customers per hour and 
for type 2 customers it is 5 customers per hour. Type 1 customers are treated with priority 
over type 2 customers. 

(i) Determine the mean sojourn time of type 1 and 2 customers under the preemptive- 
resume priority rule. 

(ii) Determine the mean sojourn time of type 1 and 2 customers under the non-preemptive 
priority rule. 

Exercise 19. 

Consider an M/M/l queue with an arrival rate of 60 customers per hour and a mean 
service time of 45 seconds. A period during which there are 5 or more customers in the 
system is called crowded, when there are less than 5 customers it is quiet. What is the 
mean number of crowded periods per day (8 hours) and how long do they last on average? 

Exercise 20. 

Consider a machine where jobs arrive according to a Poisson stream with a rate of 20 jobs 

per hour. The processing times are exponentially distributd with a mean of l//t hours. 

The processing cost is 16// dollar per hour, and the waiting cost is 20 dollar per order per 

hour. 

Determine the processing speed /* minimizing the average cost per hour. 
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Chapter 5 
M/M/c queue 



In this chapter we will analyze the model with exponential interarrival times with mean 
1/A, exponential service times with mean 1/// and c parallel identical servers. Customers 
are served in order of arrival. We suppose that the occupation rate per server, 

A 
P = — , 

C/i 

is smaller than one. 

5.1 Equilibrium probabilities 

The state of the system is completely characterized by the number of customers in the 
system. Let p n denote the equilibrium probability that there are n customers in the system. 
Similar as for the M/M/l we can derive the equilibrium equations for the probabilities p n 
from the flow diagram shown in figure 5.1. 

A A A A A 

QZEC" 603 - 

/i 2/i c/i c/i c/i 

Figure 5.1: Flow diagram for the M/M/c model 

Instead of equating the flow into and out of a single state n, we get simpler equations 
by equating the flow out of and into the set of states {0, 1, . . . , n — 1}. This amounts to 
equating the flow between the two neighboring states n — 1 and n yielding 

Xp n -i = min(n, c)/j,p n , n=l,2, ... 

Iterating gives 

(cp) n n 

Pn = — ,-Po, n = 0, ...,c 
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and 

(cp) c 
p c+n = p n Vc = p n ^~rPo, n = 0,1,2,... 
c! 

The probability po follows from normalization, yielding 

(cpT , (c P y i 



, c _! , s„ , s„ . x "I 



po= E 



\n=0 



n! c! 1 — p 



An important quantity is the probability that a job has to wait. Denote this probability 
by Iliy. It is usually referred to as the delay probability. By PASTA it follows that 

IV = p c + Pc+l + Pc+2 H 

Remark 5.1.1 (Computation ofliw) 

It will be clear that the computation of Uw by using (5.1) leads to numerical problems 
when c is large (due to the terms (cp) c /c\). In remark 11.3.2 we will formulate a numerically 
stable procedure to compute 11^^. 

5.2 Mean queue length and mean waiting time 

From the equilibrium probabilities we directly obtain for the mean queue length, 



E(L") = £>. 



c+n 



n=0 



1 r n=0 

= n w --?—, (5.2) 

1 - p 

and then from Little's law, 

E(W) = Tl w ■ . (5.3) 

1 — p Cfi 

These formulas for E(L q ) and E(W) can also be found by using the mean value technique. 
If not all servers are busy on arrival the waiting time is zero. If all servers are busy and 
there are zero or more customers waiting, then a new arriving customers first has to wait 
until the first departure and then continues to wait for as many departures as there were 
customers waiting upon arrival. An interdeparture time is the minimum of c exponential 
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(residual) service times with mean l//i, and thus it is exponential with mean 1/cp (see 
exercise 1). So we obtain 



E(W) = U w — + E(L q ) — . 
cp cp 

Together with Little's law we retrieve the formulas (5.2)-(5.3). Table 5.1 lists the delay 
probability and the mean waiting time in an M/M/c with mean service time 1 for p = 0.9. 



c 


n w 


E(W) 


1 


0.90 


9.00 


2 


0.85 


4.26 


5 


0.76 


1.53 


10 


0.67 


0.67 


20 


0.55 


0.28 



Table 5.1: Performance characteristics for the M/M/c with \x = 1 and p = 0.9 

We see that the delay probability slowly decreases as c increases. The mean waiting time 
however decreases fast (a little faster than 1/c). One can also look somewhat differently at 
the performance of the system. We do not look at the occupation rate of a machine, but at 
the average number of idle machines. Let us call this the surplus capacity. Table 5.2 shows 
for fixed surplus capacity (instead of for fixed occupation rate as in the previous table) and 
c varying from 1 to 20 the mean waiting time and the mean number of customers in the 
system. 



c 


P 


E{W) 


E(L) 


1 


0.90 


9.00 


9 


2 


0.95 


9.26 


19 


5 


0.98 


9.50 


51 


10 


0.99 


9.64 


105 


20 


0.995 


9.74 


214 



Table 5.2: Performance characteristics for the M/M/c with \x — 1 and a fixed surplus 
capacity of 0.1 server 

Although the mean number of customers in the system sharply increases, the mean 
waiting time remains nearly constant. 
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5.3 Distribution of the waiting time and the sojourn 
time 

The derivation of the distribution of the waiting time is very similar to the one in section 
4.4 for the M/M/l. By conditioning on the state seen on arrival we obtain 

oo rt+1 

P(W>t) = Y,P(J2 D k>t)Pc + n, 
n=0 fc=l 

where D^ is the fcth interdeparture time. Clearly, the random variables D^ are independent 
and exponentially distributed with mean 1/c/i. Analogously to (4.13) we find 

00 n (cut) k 

Piw>t) = EE^^V/ 

n=0 fc=0 K - 
fc=0 n=k 

p c ^ (cppt) k _ Cflt 



E~--e 



1 " P k=o kl 
= U w e- cfi{1 - p)t , t > 0. 

This yields for the conditional waiting time, 

P(W > t\W > 0) = P J^ > ^ = e - c ^- p) \ t > 0. 

Hence, the conditional waiting time W\W > is exponentially distributed with parameter 
cp(l — p). To determine the distribution of the sojourn time we condition on the length of 
the service time, so 

P(S>t) = P(W + B>t) 

roo 

= / P(W + x> t)pe~ px dx 

Jx=0 
rt roo 

= / P(W > t - x)pe~ px dx + / pe'^dx 

Jx=0 Jx=t 

= [* Il w e- c ^ l - p){t - x) pe- px dx + e- pt 

Jx=0 
= EE ( e -cv(l-p)t _ e -*\ + e -lA 

l-c(l-p) V > 



U W e -c M (l-p)t + f ! _ Uw ) e -»t 



l-c(l-p) V l-c(l-p) 

5.4 Java applet 

There is also a JAVA applet is avalaible for the performance evaluation of the M/M/c 
queue. The WWW-link to this applet is http://www.win.tue.nl/cow/Q2. 
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5.5 Exercises 

Exercise 21. (a fast and a slow machine) 

Consider two parallel machines with a common buffer where jobs arrive according to a 

Poisson stream with rate A. The processing times are exponentially distributed with mean 

1/yUi on machine 1 and l//i2 on machine 2 (/ii > ^2). Jobs are processed in order of arrival. 

A job arriving when both machines are idle is assigned to the fast machine. We assume 

that 

A 

p = — ; — 

is less than one. 

(i) Determine the distribution of the number of jobs in the system. 

(ii) Use this distribution to derive the mean number of jobs in the system. 

(iii) When is it better to not use the slower machine at all? 

(iv) Calculate for the following two cases the mean number of jobs in the system with 
and without the slow machine. 

(a) A = 2, yui = 5, fi 2 = 1; 

(b) A = 3, //1 = 5, fa = 1. 

Note: In [ ] it is shown that one should not remove the slow machine if r > 0.5 where 
r = V-il V>\- When < r < 0.5 the slow machine should be removed (and the resulting 
system is stable) whenever p < p c , where 



2 + r<1 ~ V ( 2 + r2 ) 2 + 4 (! + r2 )( 2r - !)(! + r ) 
Pc= 2(1 + r 2 ) ' 

Exercise 22. 

One is planning to build new telephone boxes near the railway station. The question is 
how many boxes are needed. Measurements showed that approximately 80 persons per 
hour want to make a phone call. The duration of a call is approximately exponentially 
distributed with mean 1 minute. How many boxes are needed such that the mean waiting 
time is less than 2 minutes? 

Exercise 23. 

An insurance company has a call center handling questions of customers. Nearly 40 calls per 
hour have to be handled. The time needed to help a customer is exponentially distributed 
with mean 3 minutes. How many operators are needed such that only 5% of the customers 
has to wait longer than 2 minutes? 

Exercise 24. 

In a dairy barn there are two water troughs (i.e. drinking places). From each trough only 
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one cow can drink at the same time. When both troughs are occupied new arriving cows 
wait patiently for their turn. It takes an exponential time to drink with mean 3 minutes. 
Cows arrive at the water troughs according to a Poisson process with rate 20 cows per 
hour. 

(i) Determine the probability that there are i cows at the water troughs (waiting or 
drinking), i = 0, 1, 2, . . . 

(ii) Determine the mean number of cows waiting at the troughs and the mean waiting 
time. 

(iii) What is the fraction of cows finding both troughs occupied on arrival? 

(iv) How many troughs are needed such that at most 10% of the cows find all troughs 
occupied on arrival? 

Exercise 25. 

A computer consists of three processors. Their main task is to execute jobs from users. 
These jobs arrive according to a Poisson process with rate 15 jobs per minute. The execu- 
tion time is exponentially distributed with mean 10 seconds. When a processor completes 
a job and there are no other jobs waiting to be executed, the processor starts to execute 
maintenance jobs. These jobs are always available and they take an exponential time with 
mean 5 seconds. But as soon as a job from a user arrives, the processor interrupts the 
execution of the maintenance job and starts to execute the new job. The execution of the 
maintenance job will be resumed later (at the point where it was interrupted). 

(i) What is the mean number of processors busy with executing jobs from users? 

(ii) How many maintenance jobs are on average completed per minute? 

(iii) What is the probability that a job from a user has to wait? 

(iv) Determine the mean waiting time of a job from a user. 
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Chapter 6 
M/E r /1 queue 



Before analyzing the M/G/l queue, we first study the M/E r /1 queue. The Erlang dis- 
tribution can be used to model service times with a low coefficient of variation (less than 
one), but it can also arise naturally. For instance, if a job has to pass, stage by stage, 
through a series of r independent production stages, where each stage takes a exponen- 
tially distributed time. The analysis of the M/E r /1 queue is similar to that of the M/M/l 
queue. 

We consider a single-server queue. Customers arrive according to a Poisson process with 
rate A and they are treated in order of arrival. The service times are Erlang-r distributed 
with mean r/jj (see subsection 2.4.4). For stability we require that the occupation rate 

p = X-- (6.1) 

is less than one. In the following section we will explain that there are two ways in which 
one can describe the state of the system. 

6.1 Two alternative state descriptions 

The natural way to describe the state of a nonempty system is by the pair (k, I) where k 
denotes the number of customers in the system and / the remaining number of service phases 
of the customer in service. Clearly this is a two-dimensional description. An alternative 
way to describe the state is by counting the total number of uncompleted phases of work 
in the system. Clearly, there is a one to one correspondence between this number and 
the pair (k,l). The number of uncompleted phases of work in the system is equal to the 
number (k — l)r + I (for the customer in service we have / phases of work instead of r). 
From here on we will work with the one- dimensional phase description. 

6.2 Equilibrium distribution 

For the one-dimensional phase description we get the flow diagram of figure 6.1. 
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Li Li il 

Figure 6.1: One-dimensional flow diagram for the M/E r /1 model 

Let p n be the equilibrium probability of n phases work in the system. By equating the 
flow out of state n and the flow into state n we obtain the following set of equilibrium 
equations for p n . 



p n {\ + Ll) 
p n (X + fi) 



Pn+iV, n = 1, 

Pn-A +Pn+lV, 



r - 1 

n = r, r + l,r + 2, ... 



(6.2) 
(6.3) 
(6.4) 



These equations may be solved as follows. We first look for solutions of (6.4) of the form 
p n = x n , n = 0,l,2,... (6.5) 



and then we construct a linear combination of these solutions also satisfying the boundary 
equations (6.2)-(6.3) and the normalization equation 



Y,Pn 
n=0 



1. 



An alternative solution approach is based on generating functions (see exercise 27). 

Substitution of (6.5) into (6.4) and then dividing by the common power x n ~ r yields the 
polynomial equation 



(A + n)x r = A + fix 



r+l 



(6.6) 



One root is x = 1, but this one is not useful, since we must be able to normalize the 
solution of the equilibrium equations. Provided condition (6.1) holds, it can be shown that 
equation (6.6) has exactly r distinct roots x with |x| < 1, say xi, . . . ,x r (see exercise 26). 
We now consider the linear combination 



Pn — z_^ CkX k , 

fc=i 



n 



0,1,2,... 



(6.7) 



For each choice of the coefficients c^ this linear combination satisfies (6.4). This freedom is 
used to also satisfy the equations (6.2)-(6.3) and the normalization equation. Note that, 
since the equilibrium equations are dependent, equation (6.2) (or one of the equations in 
(6.3)) may be omitted. Substitution of (6.7) into these equations yields a set of r linear 
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equations for r unknown coefficients. It can be shown that this set of equations has a 
unique solution, given by 

rij^fc(i - Xj/xk) 

where the subscript j runs from 1 to r. This completes the determination of the equilibrium 
probabilities p n . The important conclusion is that for the M/E r /1 queue the equilibrium 
probabilities can be expressed as a mixture of r geometric distributions. 

From the distribution of the number of phases in the system we can easily find the 
distribution of the number of customers in the system. Let qi be the probability of i 
customers in the system. Obviously qo = Po and for % > 1 it follows that 

ir 

Qi = E Pn 

n=(i — l)r J r l 

ir r 

E E c ^l 

n=(i-l)H-l fe=l 

r ir 

= E E c ^fc 

k=l n=(i-l)r+l 

k=i 

Hence, the queue length probabilities can also be expressed as a mixture of r geometric 
distributions. 

We finally remark that the roots of equation (6.6) can be numerically determined very 
efficiently (cf. [2, 3]), and that the results in this section can be easily extended to the case 
that the service times are mixtures of Erlang distributions with the same scale parameters 
(see subsection 2.4.6). 

Example 6.2.1 

Consider the M/E 2 /l queue with A = 1 and \x = 6. The equilibrium equations are given 
by 

Po = 6pi, 

7pi = 6p 2 , 

7Pn = Pn-2 + GPn+l, 71 = 2,3,4,... 

Substitution of p n = x n into the equations for n > 2 and dividing by x n ~ 2 yields 

7x 2 = 1 + 6x 3 , 
the roots of which are x = 1, x = 1/2 and x = —1/3. Hence, we set 

Pn = c 1 (-) +c 2 (--) , n = 0,1,2,... 
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and determine C\ and C2 from the equilibrium equation in n = and the normalization 
equation, 

C\ + C2 = 0C1 — 2c2, 

Cl Co 
1 = 1. 

1-1/2 1 + 1/3 

The solution is C\ = 2/5 and C2 = 4/15 (note that the equilibrium equation in n = 1 is 
also satisfied). So we obtain 

2 /l\ n 4 / l\ n 

"" = 5(2) + T5(-s) • » = °.«- 

For the queue length probabilities it follows that go = Po = 2/3 and for i > 1, 

6 /1V 8 /r ' 



?< =P2i-l+P2i = " K , . , |P . ,. 

5 \4/ 15 \9 
Note that the formula above is also valid for i = 0. 

6.3 Mean waiting time 

Let the random variable L' denote the number of phases work in the system. Then by 
PASTA it follows that 

E(W) = E(Lf)-, 
and from (6.7), 

CO 

E(L f ) = J2 n Pn 
n=0 

00 r 

= EE C K 

n=0fc=l 

r od 

= EE c k nx k 

k=l n=0 



Ck%k 

k =i (! " x k) 



Y- 



From Little's law we can also find E(L q ), the mean number of customers waiting in the 
queue. A much more direct way to determine E(W) and E(L q ) is the mean value approach. 
An arriving customer has to wait for the customers in the queue and, if the server 
is busy, for the one in service. According to the PASTA property the mean number of 
customers waiting in the queue is equal to E(L q ) and the probability that the server is 
busy on arrival is equal to p, i.e. the fraction of time the server is busy. Hence, 

E(W) = E(L q )- + pE(R), (6.8) 

A* 
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where the random variable R denotes the residual service time of the customer in service. 
If the server is busy on arrival, then with probability 1/r he is busy with the first phase 
of the service time, also with probability 1/r he is busy with the second phase, and so on. 
So the mean residual service time E(R) is equal to 



(6.9) 



E(R) 


1 r 
r /i 
r + 1 

2 


1 

■f - • 
r 

1 
1*' 


r — 1 


f ... 


1 

+ - 

r 


1 
P 


Substitution of this expression intc 


) (6.8) yields 


E(W) ■■ 


= E(L q ) L 
P 


+ P 


r + 1 
2 


1 
P ' 






Together with Little's 


law, ; 


stating 


that 






E(L q ) ■ 


= XE(W) 












we find 














E(W) ■■ 


P 
1-p 


r + 1 
2 


1 
P ' 









6.4 Distribution of the waiting time 

The waiting time can be expressed as 

Lf 

where B{ is the amount of work for the fcth phase. So the random variables Bi are inde- 
pendent and exponentially distributed with mean 1/fi. By conditioning on L* and using 
that U and Bi are independent it follows, very similar to (4.13), that 

oo n 

P{W>t) = J2 P (T, B r> t )Pn 
n=l i=l 

00 n ~ l (lltY r 

n=l i=0 l - fe=l 

fe=l i=0 n=i+l 

" fcr^fro^^ 6 

= E ^^ x fc e-^-^, t>0. 

fc=i 1 Xfc 
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Note that this distribution is a generalization of the one for the M/M/l model (namely, 
not one, but a mixture of exponentials). 

In table 6.1 we list for varying values of p and r the mean waiting time and some waiting 
time probabilities. The squared coefficient of variation of the service time is denoted by c 2 B . 
We see that the variation in the service times is important to the behavior of the system. 
Less variation in the service times leads to smaller waiting times. 



p 


r 


r 2 

( B 


E(W) 


P 


(W > 


*) 




t 5 


10 


20 


0.8 


1 


1 


4 


0.29 


0.11 


0.02 




2 


0.5 


3 


0.21 


0.05 


0.00 




4 


0.25 


2.5 


0.16 


0.03 


0.00 




10 


0.1 


2.2 


0.12 


0.02 


0.00 


0.9 


1 


1 


9 


0.55 


0.33 


0.12 




2 


0.5 


6.75 


0.46 


0.24 


0.06 




4 


0.25 


5.625 


0.41 


0.18 


0.04 




10 


0.1 


4.95 


0.36 


0.14 


0.02 



Table 6.1: Performance characteristics for the M/E r /1 with mean service time equal to 1 



6.5 Java applet 

The link to the Java applet for the performance evaluation of the M/E r /1 queue is 
http://www.win.tue.nl/cow/Q2. This applet can be used to evaluate the performance 
as a function of, e.g., the occupation rate and the shape parameter r. 
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6.6 Exercises 

Exercise 26. (Rouche's theorem) 

Rouche's theorem reads as follows. 

Let the bounded region D have as its boundary a contour C. Let the functions f(z) 

and g(z) be analytic both in D and on C, and assume that \f(z)\ < \g(z)\ on C. Then 

f(z) + g(z) has in D the same number of zeros as g(z), all zeros counted according to their 

multiplicity. 

(i) Use Rouche's theorem to prove that the polynomial equation (6.6) has exactly r 
(possibly complex) roots x with \x\ < 1. 

(Hint take f(z) = fiz r+l and g(z) = — (A + n)z r + A, and take as contour C the 
circle with center and radius 1 — e with e small and positive.) 

(ii) Show that all roots are simple. 

(Hint: Show that there are no z for which both f(z) + g(z ) and f'(z) + g'(z) vanish 
at the same time.) 

Exercise 27. (Generating function approach) 

Consider the M/E r /1 queue with arrival rate A and mean service time r//j. Let P(z) be 

the generating function of the probabilities p n , so 



'DO 



p(z) = 22 Pn z n , \z\<i. 

n=0 

(i) Show, by multiplying the equilibrium equations (6.2)-(6.4) with z n and adding over 
all states n, that 

P(z)(X + fi)- po/u = P(z)Xz r + (P(z) - p )fiz-\ 

(ii) Show that P(z) is given by 

P(z)= w {l ~ P Y r- (6-10) 

(iii) Show, by partial fraction decomposition of P(z), that the probabilities p n can be 
written as 

r / 1 \ n 

Pn = Y, C k[ — ) > n = 0,1,2,... 

where Zx, . . . , z r are the zeros of the numerator in (6.10). 
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Exercise 28. 

Orders arrive according to a Poisson process, so the inter arrival times are independent and 

exponentially distributed. Each order has to pass, stage by stage, through a series of r 

independent production stages. Each stage takes a exponentially distributed time. The 

total production time, the workload, of the order is the sum of r independent, identically 

distributed random variables. So the distribution of the workload is the r-stage Erlang 

distribution. 

(i) The state of this system (i.e., the work station together with its queue) can be 
characterized by the number of orders in the system and by the number of not yet 
completed stages for the order in the work station. 

Describe the flow diagram and give the set of equations for the equilibrium state 
probabilities. 

(ii) The state of the system at a certain time can also be described by the total number 
of production stages that have to be completed before all the orders in the system 
are ready. 

Give the flow diagram and formulate the equations for the equilibrium state proba- 
bilities. 

(iii) Define 

p n = P (total number of orders in the system = n), 

qk = P (total number of production stages in the system = k). 

Give the relation between these two probabilities. 

Exercise 29. 

Orders arrive in so called bulks. Each bulk consists of r independent orders. The bulks 

themselves arrive according to a Poisson process. The sequence in which orders of one bulk 

are processed in the work station is unimportant. All orders of a bulk have to wait until 

all orders of bulks that have arrived earlier are completed. The workload of each order is 

exponentially distributed. The state of the system is simply characterized by the number 

of orders in the system. 

Describe the flow diagram for this situation. Compare the result with that of the previous 

exercise, part (ii). 

Exercise 30. 

Jobs arrive at a machine according to a Poisson process with a rate of 16 jobs per hour. 
Each job consists of 2 tasks. Each task has an exponentially distributed processing time 
with a mean of 1 minute. Jobs are processed in order of arrival. 

(i) Determine the distribution of the number of uncompleted tasks in the system. 

(ii) Determine the distribution of the number of jobs in the system. 
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(iii) What is the mean number of jobs in the system? 
(iv) What is the mean waiting time of a job? 

Exercise 31. 

Consider an M/E 2 /l/2 queue with arrival rate 1 and mean service time 2. At time t = 

the system is empty. Determine the mean time till the first customer is rejected on arrival. 

Exercise 32. 

Consider the M/E 2 /l queue with an arrival rate of 4 customers per hour and a mean 

service time of 8 minutes. 

(i) Determine the distribution of the waiting time. 

(i) What is the fraction of customers that has to wait longer than 5 minutes? 

Exercise 33. 

Customers arrive in groups at a single-server queue. These groups arrive according to a 
Poisson process with a rate of 3 groups per hour. With probability 1/3 a group consists of 
2 customers, and with probability 2/3 it consists of 1 customer only. All customers have an 
exponential service time with a mean of 6 minutes. Groups are served in order of arrival. 
Within a group, customers are served in random order. 

(i) Determine the distribution of the number of customers in the system. 

(ii) Determine the mean sojourn time of an arbitrary customer. 

Exercise 34. 

In a repair shop of an airline company defective airplane engines are repaired. Defects occur 
according to a Poisson with a rate of 1 defective engine per 2 weeks. The mean repair time 
of an engine is 2/3 week. The repair time distribution can be well approximated by an 
Erlang-2 distribution. In the repair shop only one engine can be in repair at the same time. 

(i) Show that q n , the probability that there are n engines in the repair shop, is given by 
6 /l\ n 8 /l xn 



qn 5 \4J 15 V9 

(ii) Determine the mean sojourn time (waiting time plus repair time) of an engine in the 
repair shop. 

The airline company has several spare engines in a depot. When a defect occurs, then the 
defective engine is immediately replaced by a spare engine (if there is one available) and the 
defective engine is send to the repair shop. After repair the engine is as good as new, and 
it is transported to the depot. When a defect occurs and no spare engine is available, the 
airplane has to stay on the ground and it has to wait till a spare engine becomes available. 
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(iii) Determine the minimal number of spare engines needed such that for 99% of the 
defects there is a spare engine available. 

Exercise 35. 

Jobs arrive according to a Poisson process at a machine. The arrival rate is 25 jobs per 
hour. With probability 2/5 a job consists of 1 task, and with probability 3/5 it consists of 
2 tasks. Tasks have an exponentially distributed processing time with a mean of 1 minute. 
Jobs are processed in order of arrival. 

(i) Determine the distribution of the number of uncompleted tasks at the machine. 

(ii) Determine the mean waiting time of a job. 

(iii) Determine the mean number of jobs in the system. 

Exercise 36. 

Customers arrive in groups at a server. The group size is 2 and groups arrive according to 
a Poisson stream with a rate of 2 groups per hour. Customers are served one by one, and 
they require an exponential service time with a mean of 5 minutes. 

(i) Determine the distribution of the number of customers in the system. 

(ii) What is the mean waiting time of the first customer in a group? 

(iii) And what is the mean waiting time of the second one? 
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Chapter 7 

M/G/l queue 



In the M/G/l queue customers arrive according to a Poisson process with rate A and 
they are treated in order of arrival. The service times are independent and identically 
distributed with distribution function Fb(-) and density /b(-). For stability we have to 
require that the occupation rate 

p = XE(B) (7.1) 

is less than one. In this chapter we will derive the limiting or equilibrium distribution 
of the number of customers in the system and the distributions of the sojourn time, the 
waiting time and the busy period duration. It is further shown how the means of these 
quatities can be obtained by using the mean value approach. 

7.1 Which limiting distribution? 

The state of the M/G/l queue can be described by the pair (n,x) where n denotes the 
number of customers in the system and x the service time already received by the customer 
in service. We thus need a two-dimensional state description. The first dimension is still 
discrete, but the other one is continuous and this essentially complicates the analysis. 
However, if we look at the system just after departures, then the state description can be 
simplified to n only, because x = for the new customer (if any) in service. Denote by Lf 
the number of customers left behind by the fcth departing customer. In the next section 
we will determine the limiting distribution 

d n = lim P(L d k = n). 

The probability d n can be interpreted as the fraction of customers that leaves behind n 
customers. But in fact we are more interested in the limiting distribution p n defined as 

p n = lim P(L(t) = n), 

t — >oo 

where L(t) is the number of customers in the system at time t. The probability p n can 
be interpreted as the fraction of time there are n customers in the system. From this 
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distribution we can compute, e.g., the mean number of customers in the system. Another 
perhaps even more important distribution is the limiting distribution of the number of 
customers in the system seen by an arriving customer, i.e., 

a n = lim P(L a k = n), 

k— >oo 

where L k is the number of customers in the system just before the fcth arriving customer. 
From this distribution we can compute, e.g., the distribution of the sojourn time. What is 
the relation between these three distributions? It appears that they all are the same. 

Of course, from the PASTA property we already know that a n = p n for all n. We will 
now explain why also a n = d n for all n. Taking the state of the system as the number of 
customers therein, the changes in state are of a nearest-neighbour type: if the system is in 
state n, then an arrival leads to a transition from n to n + 1 and a departure from n to 
n—1. Hence, in equilibrium, the number of transitions per unit time from state n to n + 1 
will be the same as the number of transitions per unit time from n + 1 to n. The former 
transitions correspond to arrivals finding n customers in the system, the frequency of which 
is equal to the total number of arrivals per unit time, A, multiplied with the fraction of 
customers finding n customers in the system, a n . The latter transitions correspond to 
departures leaving behind n customers. The frequency of these transitions is equal to the 
total number of departures per unit time, A, multiplied with the fraction of customers 
leaving behind n customers, d n . Equating both frequencies yields a n = d n . Note that 
this equality is valid for any system where customers arrive and leave one by one (see also 
exercise 48). Thus it also holds for, e.g., the G/G/c queue. 

Summarizing, for the M/G/l queue, arrivals, departures and outside observers all see 
the same distribution of number of customers in the system, i.e., for all n, 

7.2 Departure distribution 

In this section we will determine the distribution of the number of customers left behind 
by a departing customer when the system is in equilibrium. 

Denote by Lf the number of customers left behind by the kth departing customer. We 
first derive an equation relating the random variable L k+l to Lf. The number of customers 
left behind by the k + 1th customer is clearly equal to the number of customers present 
when the kth. customer departed minus one (since the k + 1th customer departs himself) 
plus the number of customers that arrives during his service time. This last number is 
denoted by the random variable A k +i- Thus we have 

-kfc+i = L k - 1 + A k+ x, 
which is valid if Lf > 0. In the special case Lf = 0, it is readily seen that 

L k+1 = A k+1 . 
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50 forms a 



From the two equations above it is immediately clear that the sequence {Lf}% 
Markov chain. This Markov chain is usually called the imbedded Markov chain, since we 
look at imbedded points on the time axis, i.e., at departure instants. 
We now specify the transition probabilities 



P 



i j 



P(L d k+1 =j\L d k = i). 

Clearly pij = for all j < i — 1 and pij for j > i — 1 gives the probability that exactly 
j — % + 1 customers arrived during the service time of the k + 1th customer. This holds 
for i > 0. In state the kth customer leaves behind an empty system and then p j gives 
the probability that during the service time of the k + 1th customer exactly j customers 
arrived. Hence the matrix P of transition probabilities takes the form 

\ 



P 



J 

where a n denotes the probability that during a service time exactly n customers arrive. 
To calculate a n we note that given the duration of the service time, t say, the number of 
customers that arrive during this service time is Poisson distributed with parameter At. 
Hence, we have 

f°° (Xt) n 
a n = / — j- 
Jt=o n\ 

The transition probability diagram is shown in figure 7.1. 
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Figure 7.1: Transition probability diagram for the M/G/l imbedded Markov chain 

This completes the specification of the imbedded Markov chain. We now wish to 
determine its limiting distribution. Denote the limiting distribution of Lf. by {d n }%L and 
the limiting random variable by L d . So 



d n = P(L° 



n 



lim P{L% 

k— >oo 



n\ 
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The limiting probabilities d n , which we know are equal to p n , satisfy the equilibrium 
equations 

d n = rf„+i«o + d n a,i + • • • + d\a n + d a n 

n 

= Y dn+i-kUk + d a n , n = 0, 1, . . . (7.3) 

fc=0 

To solve the equilibrium equations we will use the generating function approach. Let us 
introduce the probability generating functions 

oo oo 

P L a{ Z ) = Y d n Z n , P A (Z) = E a nZ n , 

n=0 n=0 

which are defined for all z < 1. Multiplying (7.3) by z n and summing over all n leads to 

oo / n \ 

Pl4 z ) = Yl ( Y d n+ i- k a k + d a n J z n 

n=0 \fc=0 / 

oo n oo 

= z"'EE d n+1 . k z n+1 - k a k z k + Y do® n z n 

n=0 k=0 n=0 

oo oo 

= z- l YY. d n+i-kZ n+1 - k a k z k + d Q P A {z) 

k=0 n=k 

oo oo 

= Z- l Y*kZ k Y d n + l-kZ n+l - k + d P A ( Z ) 
A;=0 n=k 

= z- 1 P A (z)(P Ld (z) - d ) + d P A {z). 

Hence we find 

d P A (z)(l - z- 1 ) 



P Ld {z) 



1 - z-*P A {z) 

To determine the probability do we note that do is equal to po, which is the fraction of time 
the system is empty. Hence do = po = 1 — p ( alternatively, do follows from the requirement 
P L d(l) = 1). So, by multiplying numerator and denominator by —z we obtain 

By using (7.2), the generating function P A (z) can be rewritten as 

P A {z) = Y [ -^-e- xt f B (t)dtz n 

^ Jt=o n\ 

■I i=o ^ n\ 

e- {x - Xz)t f B (t)dt 

t=Q 

B{X-\z) (7.5) 
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Substitution of (7.5) into (7.4) finally yields 

L k ' B(A - Az) - z * ' 

This formula is one form of the Pollaczek-Khinchin formula. In the following sections 
we will derive similar formulas for the sojourn time and waiting time. By differentiating 
formula (7.6) we can determine the moments of the queue length (see section 2.2). To find 
its distribution, however, we have to invert formula (7.6), which usually is very difficult. In 
the special case that B(s) is a quotient of polynomials in s, i.e., a rational function, then in 
principle the right-hand side of (7.6) can be decomposed into partial fractions, the inverse 
transform of which can be easily determined. The service time has a rational transform 
for, e.g., mixtures of Erlang distributions or Hyperexponential distributions (see section 
2.4). The inversion of (7.6) is demonstrated below for exponential and Erlang-2 service 
times. 

Example 7.2.1 (M/M/l) 

Suppose the service time is exponentially distributed with mean 1/p. Then 

B W --e-. 

p + s 
Thus 

_ (i-p) M rfa( 1 - 2 ) _ (1-pMi-z) _ (i-p) f i(i-z) _ l-p 

Ld[Z) j+t^-z p-z(p + X-Xz) (p-Xz)(l-z) 1-pz' 

Hence 

d n = Vn = (1 - p)p n , n = 0,1,2,... 

Example 7.2.2 (M/E 2 /l) 

Suppose the service time is Erlang-2 distributed with mean 2/ p. Then (see subsection 
2.4.4) 



so 



B(s) 



P Ld (z) 




2 
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P) V^+A-Aj 


-Ya-z) 
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\/a+X-XzJ 

--p)p 2 (l- 


-z) 


p?- 


- z(p + A - 

(l-p)(l- 


Xz) 2 

-z) 


1- 


z(l + p(l- 
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-z)/2)* 



1 - pz - p 2 z(l - z)/4 ' 
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For p = 1/3 we then find 

2/3 24 



Pia{z) 



l-z/3-z(l-z)/36 36-13z + z 2 

24 24/5 24/5 6/5 8/15 



(4-z)(9-z) 4-2 9-2 l-z/4 1 — jzr/9 

Hence, 

d B =p B = -( i J --(-) , n = 0,l,2,... 

Example 7.2.3 (M/H 2 /l) 

Suppose that A = 1 and that the service time is hyperexponentially distributed with 
parameters pi = 1 — P2 = 1/4 and pi = 1, p 2 = 2. So the mean service time is equal to 
1/4 • 1 + 3/4 • 1/2 = 5/8. The Laplace- Stieltjes transform of the service time is given by 
(see subsection 2.4.5) 

~, N 1 1 3 2 1 8 + 7s 

Bis) = - • 1 = - • ~, -, r . 

y ' 4 1 + s 4 2 + s 4 (l + s)(2 + s) 
Thus we have 

3 1 15-7, _ {l _ z) 



P Ld (z) 



8 4 (2-z)(3-z) ' 



1 15-7, 



4 (2-z)(3-*) 

3 (lb-7z)(l-z) 

8 ' {15-7z)-Az(2-z)(3-z) 

3 15-7z 



8 (3-2z)(5-2z) 
3 9/4 3 5/4 



8 3-22 8 5-2^/5 

9/32 3/32 

1 - 22/3 + 1 - 22/5 ' 



So 



9 (2\ n 3 (2\ n 
dn=Pn = ^{- 3 ) +^y , n = 0,l,2,... 

7.3 Distribution of the sojourn time 

We now turn to the calculation of how long a customer spends in the system. We will show 
that there is a nice relationship between the transforms of the time spent in the system 
and the departure distribution. 

Let us consider a customer arriving at the system in equilibrium. Denote the total 
time spent in the system for this customer by the random variable S with distribution 
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function Fs(-) and density /s(-). The distribution of the number of customers left behind 
upon departure of our customer is equal to {d n }^ =0 (since the system is in equilibrium). 
In considering a first-come first-served system it is clear that all customers left behind are 
precisely those who arrived during his stay in the system. Thus we have (cf. (7.2)) 

(Ai)\_ Ai 



d n = / ^-e- M f s (t)dt. 
Jt=o n! 

Hence, we find similarly to (7.5) that 

p Ld ( z ) = S(X-Xz). 
Substitution of this relation into (7.6) yields 

B{\-\z)-z 
Making the change of variable s = A — Xz we finally arrive at 

XB{s) + s — A 
This formula is also known as the Pollaczek-Khinchin formula. 

Example 7.3.1 (M/M/l) 

For exponential service times with mean \j\x we have 

B(») = -*-. 

Thus 



S(s) 



(i-p)^s __ (i-pVs (i-pVs Mi-p) 



A^ + s-A A/i + (s- A)(/i + s) (/i-A)s + s 2 yu(l-p) + s' 

Hence, 5 is exponentially distributed with parameter nil — p), i.e., 

F s (t) = P(S <t) = l- e"^ 1 -^, t > 0. 

Example 7.3.2 (M/E 2 /l) 

Suppose that A = 1 and that the service time is Erlang-2 distributed with mean 1/3, so 

sw ' 6 2 



,6 + s 
Then it follows that (verify) 

F s (t) = 8 5 (l-e- 3t )- 3 5 (l-e- st ), t>0. 

Example 7.3.3 (M/H 2 /l) 

Consider example 7.2.3 again. From (7.7) we obtain (verify) 

F s (t) = — (1 - e -'/2) + 1 (1 _ e -3t/2) t > Q 
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7.4 Distribution of the waiting time 

We have that S, the time spent in the system by a customer, is the sum of W (his waiting 
time) and B (his service time), where W and B are independent. Since the transform of 
the sum of two independent random variables is the product of the transforms of these two 
random variables (see section 2.3), it holds that 

S(s) = W{s)-B(s). (7.8) 

Together with (7.7) it follows that 

W(s) = J-f , (7.9) 

XB(s) + s — A 

which is the third form of the Pollaczek-Khinchin formula. 

Example 7.4.1 (M/M/l) 

For exponential service times with mean 1/p we have 

B(s) " 



P + S 

Then from (7.9) it follows that (verify) 

fi{l-p) + s 

The inverse transform yields 

F w (t) = P{W < t) = (1 - p) + p(l - e-^- p)t ), t > 0. 

Hence, with probability (1 — p) the waiting time is zero (i.e., the system is empty on arrival) 
and, given that the waiting time is positive (i.e., the system is not empty on arrival), the 
waiting time is exponentially distributed with parameter p(l — p) (see also section 4.4). 

7.5 Lindley's equation 

In the literature there are several alternative approaches to derive the Pollaczek-Khinchin 
formulas. In this section we present one that is based on a fundamental relationship 
between the waiting time of the nth customer and the n + 1th customer. 

Denote by A n the interarrival time between the nth and n + lth customer arriving at the 
system. Further let W n and B n denote the waiting time and the service time, respectively, 
for the nth customer. Clearly, if A n < W n + B n , then the waiting time of the n + 1th 
customer is equal to W n + B n — A n . If on the other hand A n > W n + B n , then his waiting 
time is equal to zero. Summarizing, we have 

W n+1 = max(H/ n + B n - A n , 0). (7.10) 
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This equation is commonly referred to as Lindley's equation. We denote the limit of W n 
as n goes to infinity by W (which exists if p is less than one). Then, by letting n go to 
infinity in (7.10) we obtain the limiting form 

W = mnx(W + B-A,0) (7.11) 

(where we also dropped the subscripts of A n and B n , since we are considering the limiting 
situation). With S = W + B, this equation may also be written as 

W = max(S - A, 0) (7.12) 

Note that the interarrival time A is exponentially distributed with parameter A and that 
the random variables S and A are independent. We now calculate the transform of W 
directly from equation (7.12). This gives 

W(s) = E(e- sW ) 

_ P(e~ smax (S— A,0)\ 
oo roo 

/ e- smax{x - y ' 0) f s (x)\e- Xy dxdy 

x=Q Jy=0 

OO PQj P OO /*00 

/ e- s{x - y) f s {x)\e- Xy dxdy + / / f s {x)\e- Xy dxdy 

x=0 Jy=0 J x=0 Jy=x 

\ roo roo 

— Ax „— sx\ r / \ i i / „— \x 



( e -Ax _ e s*)f s ( x ) dx + / e- Ax f s {x)dx 



S — X Jx=0 Jx=0 

= -^(S(X)-S(s)) + S(X) 
s — X 

Substitution of (7.8) then yields 

XB(s) + s — X 
It remains to find S(X). Realizing that 

S(X) = / e- Xx f s (x)dx = / P(A > x)f s (x)dx = P{A > S), 

it follows that S(X) is precisely the probability that the system is empty on arrival, so 

S(X) = 1- p. 
Substitution of this relation into (7.13) finally yields formula (7.9). 
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7.6 Mean value approach 

The mean waiting time can of course be calculated from the Laplace-Stieltjes transform 
(7.9) by differentiating and substituting s = (see section 2.3). In this section we show 
that the mean waiting time can also be determined directly (i.e., without transforms) with 
the mean value approach. 

A new arriving customer first has to wait for the residual service time of the customer 
in service (if there is one) and then continues to wait for the servicing of all customers who 
were already waiting in the queue on arrival. By PASTA we know that with probability p 
the server is busy on arrival. Let the random variable R denote the residual service time 
and let L q denote the number of customers waiting in the queue. Hence, 

E(W) = E(L q )E(B) + pE(R), 
and by Little's law (applied to the queue), 

E(L q ) = XE(W). 
So we find 

E{W) = pE ^. (7.14) 

1- p 

Formula (7.14) is commonly referred to as the Pollaczek-Khinchin mean value formula. It 
remains to calculate the mean residual service time. In the following section we will show 
that 

which may also be written in the form 

E ^-- E ^r al ^--\ {A+l)E{B) - (7 - 16) 

An important observation is that, clearly, the mean waiting time only depends upon the 
first two moments of service time (and not upon its distribution). So in practice it is 
sufficient to know the mean and standard deviation of the service time in order to estimate 
the mean waiting time. 

7.7 Residual service time 

Suppose that our customer arrives when the server is busy and denote the total service 
time of the customer in service by X . Further let fx(') denote the density of X. The basic 
observation to find fx(') is that it is more likely that our customer arrives in a long service 
time than in a short one. So the probability that X is of length x should be proportional 
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to the length x as well as the frequency of such service times, which is fB(x)dx. Thus we 
may write 

P(x < X < x + dx) = fx(x)dx = Cxfs(x)dx, 

where C is a constant to normalize this density So 

poo 

C~ l = / xf B {x)dx = E(B). 

Hence 

, , v xf B (x) 

Given that our customer arrives in a service time of length x, the arrival instant will be 
a random point within this service time, i.e., it will be uniformly distributed within the 
service time interval (0,x). So 

dt 
P(t< R<t + dt\X = x) = — , t<x. 

x 

Of course, this conditional probability is zero when t > x. Thus we have 

r°° dt r°° fa(r) 1 — Fa(t) 

P(t<R<t + dt) = f R (t)dt = / -fx(x)dx = / ^U. dxdt = * U dt. 

Jx=t X Jx=t -Cj{-D) ^{^) 

This gives the final result 

f (t) - 1 -— ^ 
fR{t > ~ E(B) ' 

from which we immediately obtain, by partial integration, 

E{B 2 ) 



roo J roo I roo J 

b(h) = / t M m = scb) to f(1 - F »w>* = Em Lf ut)dt 



E(B)Jt=o x w/ E{B)Jt=o2 JDW 2E(B) 

This computation can be repeated to obtain all moments of R, yielding 
E{B n+1 ) 



E(R n ) 



(n+l)E(B) 



Example 7.7.1 (Erlang service times) 

For an Erlang-r service time with mean r//j we have 



E(B) = -, o\B) = - 2 



i r :» ;7 , l , tT?tU ^-i = r{l+r) 

2 



SO 

E(B 2 )=a 2 (B)r(E(B)) 

Hence (see also (6.9)) 
1 + r 

E ^ - ir 
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7.8 Variance of the waiting time 

The mean value approach fails to give more information than the mean of the waiting time 
(or other performance characteristics). But to obtain information on the variance or higher 
moments of the waiting time we can use formula (7.9) for the Laplace-Stieltjes transform 
of the waiting time. This is demonstrated below. 
We first rewrite formula (7.9) as 
1-p 



W(s) 



1 P I sE{B) 



The term between brackets can be recognised as the transform of R, since by partial 
integration we have 

R(s) = E(e- sR ) = / e- st f R (t)dt = — — / e~ st (l - F B (t))dt 

Jt=0 E{B) Jt=0 



EiB){-s-Ls e fBm ) = ^E{Bj 



Hence 



w(s) - Y^k ■ (7 - 17) 

By differentiating (7.17) and substituting s = we retrieve formula (7.14), i.e., 

1-p 

By differentiating (7.17) twice and substituting s = we find 

E(W 2 ) = 2(E(W)) 2 + pE ^. 

1-p 

Hence, for the variance of the waiting time we now obtain 

a\W) = E(W 2 ) - {E{W)f = (I ^ o)2 {p{E{R)f + (1 - p)E(R 2 )) . 

The first two moments of the conditional waiting time W\W > can be calculated 
from 

E(W\W > 0) = ^^^ , E(W 2 \W > 0) = E ^ W ^ . 
P P 

It then follows, after some algebra, that the squared coefficient of variation of W|W > 

is given by the simple equation 

c w\w>o = p + (1 - p)4 = 1 + (i - p)(4 - 1). 

This implies that for p near to 1 the squared coefficient of variation of the conditional 
waiting time is close to 1. Hence, as a rule of thumb, we may think of the conditional 
waiting time as an exponential random variable (which has a coefficient of variation equal 
to 1) and thus obtain a rough estimate for its distribution. 
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7.9 Distribution of the busy period 

The mean duration of a busy period in the M/G/l queue can be determined in exactly 
the same way as for the M/M/l queue in subsection 4.6.1. Thus we have 

E(BP) = 1 ^. (7.18) 

1- p 

The calculation of the distribution of a busy period in the M/G/l queue is more compli- 
cated. 

We first derive an important relation for the duration of a busy period. What happens 
in a busy period? We start with the service of the first customer. His service time is 
denoted by B\. During his service new customers may arrive to the system. Denote this 
(random) number by TV and label these customers by Ci,...,Cjv. At the departure of 
the first customer we take customer C\ into service. However, instead of letting the other 
customers C2, . . . , Cn wait for their turn, we choose to take them temporarily out of the 
system. Customer C2 will be put into the system again and taken into service as soon as 
the system is empty again. So customers arriving during the service of C\ will be served 
first (as well as the ones arriving during their service, and so on). Thus it is as if C\ initiates 
a new busy period for which C2 has to wait. This busy period will be called a sub-busy 
period. In the same way C3 has to wait for the sub-busy period initiated by C2, and so 
on. Finally the sub-busy period due to Cn completes the major busy period. Note that 
our customers C2, . . . , Cn are not treated first-come first-served anymore. But this does 
not affect the duration of the (major) busy period, because its duration is independent of 
the order in which customers are served. Denote the busy period by BP and the sub-busy 
period due to Cj by BPi. Then we have the important relation 

BP = B x + BP 1 + • • • + BP N (7.19) 

where the random variables BP, BP\, BP2, . . . are independent and all have the same dis- 
tribution, and further, they are independent of N. The number TV of course depends on 
the service time B\. 

We will now use relation (7.19) to derive the Laplace-Stieltjes transform of the busy 
period BP. By conditioning on the length of B\ we get 

roc 

BP{s) = E{e~ sBP ) = / E{e~ sBP \B l = t)f B {t)dt 

Jt=o 

and by also conditioning on N, 

00 
E{e~ sBP \B l = t) = Y, E ( e ~ sBP \ B ± = t ' N = n)P(N = n\B x = t) 

n=0 

= f; E(e-& +BPl+ - +Bp »)\B 1 =t,N = n)^-e 



(Ai)\_ Ai 



n=0 
00 



ill 



(t+BP 1 +-+BP n )\ \™r_ -\t 



J2 E ^sit+BP 1+ - + BP n )- ) __ ( 



n=0 
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f; E{e- St ■ e- sBP ' ■ ■ ■ e sBP n) (^L e -xt 

n=0 n - 



OG 



^ e-°\E(e-° BP )) 
=o 

-{s+\-\BP(s))t 



,.IJl>,,n\"t) -\t 



n=0 U - 



Thus we have 

BP{s) = / e^^^V^R 

Jt=0 

which can be recognised as 

BP(s) = B(s + X-XBP(s)). (7.20) 

In the example below it is shown that even for the simple M/M/l queue it is not easy to 
determine the distribution of BP from this equation. 

Example 7.9.1 (M/M/l) 

For the M/M/l we have 

B(s) = -S- . 

fl + S 

In this case equation (7.20) reduces to 

BP(s) = ^ == . 

fi + s + X- XBP(s) 

So BP(s) is a root of the quadratic equation 

X(BP(s)) 2 - (A + fi + s)BP(s) + ix = 0. 

Solving this equation and restricting our solution to the case for which < BP(s) < 1 for 
all s > yields 

1 



£P(s) = 2X ( A + ^ + s ~ V ( A + ^ + s ) 2 " 4A ^J • 

This equation may be inverted to obtain the density of the busy period. This gives an 
expression involving a Bessel function (see subsection 4.6.2). 

It is, however, straightforward to determine the moments of BP from equation (7.20) 
by differentiating and substituting s = (see section 2.3). This will be demonstrated 
below. 

Differentiating (7.20) and substituting s = yields 

-E(BP) = BP W (0) = B {1) (0) ■ (1 - XBP (1 \o)) = -E(B)(1 + XE(BP)). 
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Hence, with p = \E(B), we retrieve formula (7.18), i.e., 

E(BP) = ' E ^-. (7.21) 

1- p 

By differentiating (7.20) twice and substituting s = we get 

E(BP 2 ) = BP {2 \0) 

= 5 (2) (0) • (1 - XBP (1 \o)) 2 + 5 (1) (0) • (-XBP {2 \0)) 
= E(B 2 )(1 + XE(BP)) 2 + XE(B)E(BP 2 ). 

From this we obtain 

E(BP 2 ) = 0^ . (7.22) 

From (7.21) and (7.22) it follows that the squared coefficient of variation, c 2 BP , is given by 



4 , p 



^-i-p^i- p - 

Clearly, as p tends to one, the variation in BP explodes. 



Example 7.9.2 (M/M/l) 

Since c 2 B = 1 for exponential service times it follows from (7.23) that 

r 2 _i±P 
CBP ~ 1 - P • 

In table 7.1 we list c 2 BP for increasing values of p. 



(7.23) 



p 


0.5 


0.6 


0.7 


0.8 


0.9 


0.95 


r 2 
C BP 


3 


4 


5.7 


9 


19 


39 



Table 7.1: The squared coefficient of the busy period for the M/M/l queue 



7.10 Java applet 

There is a JAVA applet avalaible for the evaluation of mean performance characteristics 
of the M/G/l queue. The link to this applet is http://www.win.tue.nl/cow/Q2. 
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7.11 Exercises 

Exercise 37. (Post office) 

At a post office customers arrive according to a Poisson process with a rate of 30 customers 
per hour. A quarter of the customers wants to cash a cheque. Their service time is 
exponentially distributed with a mean of 2 minutes. The other customers want to buy 
stamps and their service times are exponentially distributed with a mean of 1 minute. In 
the post office there is only one server. 

(i) Determine the generating function Pl(z) of the number of customers in the system. 

(ii) Determine the distribution of the number of customers in the system. 
(iii) Determine the mean number of customers in the system, 
(iv) Determine S(s), the Laplace-Stieltjes transform of the sojourn time. 

(v) Determine the mean and the distribution of the sojourn time. 

(vi) Determine the mean busy period duration. 

(vii) Determine the mean number of customers in the system and the mean sojourn time 
in case all customers have an exponentially distributed service time with a mean of 
75 seconds. 

Exercise 38. 

Consider a single machine where jobs arrive according to a Poisson stream with a rate of 
10 jobs per hour. The processing time of a job consists of two phases. Each phase takes 
an exponential time with a mean of 1 minute. 

(i) Determine the Lapace-Stieltjes transform of the processing time. 

(ii) Determine the distribution of the number of jobs in the system. 

(iii) Determine the mean number of jobs in the system and the mean production lead 
time (waiting time plus processing time). 

Exercise 39. (Post office) 

At a post office customers arrive according to a Poisson process with a rate of 60 customers 

per hour. Half of the customers have a service time that is the sum of a fixed time of 15 

seconds and an exponentially distributed time with a mean of 15 seconds. The other half 

have an exponentially distributed service time with a mean of 1 minute. 

Determine the mean waiting time and the mean number of customers waiting in the queue. 

Exercise 40. (Two phase production) 

A machine produces products in two phases. The first phase is standard and the same for 

all products. The second phase is customer specific (the finishing touch). The first (resp. 
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second) phase takes an exponential time with a mean of 10 (resp. 2) minutes. Orders for 
the production of one product arrive according to a Poisson stream with a rate of 3 orders 
per hour. Orders are processed in order of arrival. 
Determine the mean production lead time of an order. 

Exercise 41. 

Consider a machine where jobs are being processed. The mean production time is 4 minutes 

and the standard deviation is 3 minutes. The mean number of jobs arriving per hour is 

10. Suppose that the inter arrival times are exponentially distributed. 

Determine the mean waiting time of the jobs. 

Exercise 42. 

Consider an M/G/l queue, where the server successfully completes a service time with 
probability p. If a service time is not completed successfully, it has to be repeated until it 
is successful. Determine the mean sojourn time of a customer in the following two cases: 

(i) The repeated service times are identical. 
(ii) The repeated service times are independent, and thus (possibly) different. 

Exercise 43. 

At the end of a production process an operator manually performs a quality check. The 

time to check a product takes on average 2 minutes with a standard deviation of 1 minute. 

Products arrive according to a Poisson stream. One considers to buy a machine that is able 

to automatically check the products. The machine needs exactly 84 seconds to perform a 

quality check. Since this machine is expensive, one decides to buy the machine only if it 

is able to reduce lead time for a quality check (i.e. the time that elapses from the arrival 

of a product till the completion of its quality check) to one third of the lead time in the 

present situation. 

So only if the arrival rate of products exceeds a certain threshold, one will decide to buy 

the machine. Calculate the value of this threshold. 

Exercise 44. (Robotic dairy barn) 

In a robotic dairy barn cows are automatically milked by a robot. The cows are lured into 
the robot by a feeder with nice food that can only be reached by first passing through 
the robot. When a cow is in the robot, the robot first detects whether the cow has to 
be milked. If so, then the cow will be milked, and otherwise, the cow can immediately 
leave the robot and walk to the feeder. A visit to the robot with (resp. without) milking 
takes an exponential time with a mean of 6 (resp. 3) minutes. Cows arrive at the robot 
according to a Poisson process with a rate of 10 cows per hour, a quarter of which will be 
milked. 

(i) Show that the Laplace-Stieltjes transform of the service time in minutes of an arbi- 
trary cow at the robot is given by 

B(s) ' 4 + 21S 



4 (l + 6s)(l + 3s) 
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(ii) Show that the Laplace-Stieltjes transform of the waiting time in minutes of an arbi- 
trary cow in front of the robot is given by 

vTr f \ 3 9 l 1 l 

W(s) = - 



8 16 l + 12s 16 l + 4s 

(iii) Determine the fraction of cows for which the waiting time in front of the robot is less 
than 3 minutes. 

(iv) Determine the mean waiting time in front of the robot. 

Exercise 45. 

Consider a machine processing parts. These parts arrive according to a Poisson stream 
with a rate of 1 product per hour. The machine processes one part at a time. Each part 
receives two operations. The first operation takes an exponential time with a mean of 15 
minutes. The second operation is done immediately after the first one and it takes an 
exponential time with a mean of 20 minutes. 

(i) Show that the Laplace-Stieltjes transform of the production lead time (waiting time 
plus processing time) in hours is given by 

-51 15 

S(s) 



4 1+s 4 5+s 

(ii) Determine the distribution of the production lead time and the mean production lead 
time. 

One uses 3 hours as a norm for the production lead time. When the production lead time 
of a part exceeds this norm, it costs 100 dollar. 

(iii) Calculate the mean cost per hour. 

Exercise 46. (Warehouse) 

In a warehouse for small items orders arrive according to a Poisson stream with a rate of 6 
orders per hour. An order is a list with the quantities of products requested by a customer. 
The orders are picked one at a time by one order picker. For a quarter of the orders the 
pick time is exponentially distributed with a mean of 10 minutes and for the other orders 
the pick time is exponentially distributed with a mean of 5 minutes. 

(i) Show that the Laplace-Stieltjes transform of the pick time in minutes of an arbitrary 
order is given by 

g(.) x 4+35s 



4 (l + 10s)(l + 5s) 
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(ii) Show that the Laplace-Stieltjes transform of the lead time (waiting time plus pick 
time) in minutes of an arbitrary order is given by 

~ 5 3 27 1 

S(s) 



32 3 + 20s 32 1 + 20s 
(iii) Determine the fraction of orders for which the lead time is longer than half an hour. 
(iv) Determine the mean lead time. 

Exercise 47. (Machine with breakdowns) 

Consider a machine processing parts. Per hour arrive according to a Poisson process on 
average 5 orders for the production of a part. The processing time of a part is exactly 
10 minutes. During processing, however, tools can break down. When this occurs, the 
machine immediately stops, broken tools are replaced by new ones and then the machine 
resumes production again. Replacing tools takes exactly 2 minutes. The time that elapses 
between two breakdowns is exponentially distributed with a mean of 20 minutes. Hence, 
the total production time of a part, B, consists of the processing time of 10 minutes plus 
a random number of interruptions of 2 minutes to replace broken tools. 

(i) Determine the mean and variance of the production time B. 
(ii) Determine the mean lead time (waiting time plus production time) of an order. 

Exercise 48. (Arrival and departure distribution) 

Consider a queueing system in which customers arrive one by one and leave one by one. 
So the number of customers in the system only changes by +1 or — 1. We wish to prove 
that a n = d n . Suppose that the system is empty at time t = 0. 

(i) Show that if Lf +1 < n, then L% +n+1 < n. 
(ii) Show that if Ll +n+1 < n, then L d k < n. 
(iii) Show that (i) and (ii) give, for any n > 0, 

lim P(L d k <n)= lim P(L a k < n). 

fc^oo k— >oo 

Exercise 49. (Number served in a busy period) 

Let the random variable Nb p denote the number of customers served in a busy period. 

Define the probabilities /„ by 

f n = P(N bp = n). 

We now wish to find its generating function 

oo 

F Nbp (z) = J2f n z n . 



n=l 
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(i) Show that F^ b (z) satisfies 

F Nbp (z) = zB(X-XF Nbp (z)). (7.24) 

(ii) Show that the mean number of customers served in a busy period is given by 

E(N bp ) = ^i- . 

(iii) Solve equation (7.24) for the M/M/l queue. 

(iv) Solve equation (7.24) for the M/D/l queue (Hint: Use Lagrange inversion formula, 
see, e.g., [30]). 
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Chapter 8 

G/M/l queue 



In this chapter we study the G/M/l queue, which forms the dual of the M/G/l queue. In 
this system customers arrive one by one with interarrival times identically and indepen- 
dently distributed according to an arbitrary distribution function Fa(-) with density /a(0- 
The mean interarrival time is equal to 1/A. The service times are exponentially distributed 
with mean l//x. For stability we again require that the occupation rate p = X/fi is less 
than one. 

The state of the G/M/l queue can be described by the pair (n, x) where n denotes the 
number of customers in the system and x the elapsed time since the last arrival. So we 
need a complicated two-dimensional state description. However, like for the M/G/l queue, 
the state description is much easier at special points in time. If we look at the system on 
arrival instants, then the state description can be simplified to n only, because x = at an 
arrival. Denote by L% the number of customers in the system just before the fcth arriving 
customer. In the next section we will determine the limiting distribution 

a n = lim P{L% = n). 

k— >oo 

From this distribution we will be able to calculate the distribution of the sojourn time. 

8.1 Arrival distribution 

In this section we will determine the distribution of the number of customers found in the 
system just before an arriving customer when the system is in equilibrium. 

We first derive a relation between the random variables L k+1 and L k . Defining the 
random variable D k+ i as the number of customers served between the arrival of the kth 
and k + 1th customer, it follows that 

L a k+l = L a k + 1 - D k+1 . 

From this equation it is immediately clear that the sequence {Lfy^Q forms a Markov chain. 
This Markov chain is called the G/M/l imbedded Markov chain. 
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We must now calculate the associated transition probabilities 



Clearly pij = for all j > i + 1 and pij for j < i-\- 1 is equal to the probability that exactly 
i + 1 — j customers are served during the interarrival time of the k + 1th customer. Hence 
the matrix P of transition probabilities takes the form 

\ 



P 



' Po,o 


00 









Pl,0 


0i 


00 







P2,0 


02 


01 


00 





P3,0 


03 


fh 


01 


00 



V 



/ 



where 0i denotes the probability of serving % customers during an interarrival time given 
that the server remains busy during this interval (thus there are more than % customers 
present). To calculate 0i we note that given the duration of the interarrival time, t say, 
the number of customers served during this interval is Poisson distributed with parameter 
fit. Hence, we have 



0i 



(^) „-/rf 



e-^f A (t)dt. 



h=o i\ 
Since the transition probabilities from state j should add up to one, it follows that 



(8.1) 



P^fi = l-£&= J2 03- 

3=0 j=i+l 

The transition probability diagram is shown in figure 8.1. 



A 




Figure 8.1: Transition probability diagram for the G/M/l imbedded Markov chain 

This completes the specification of the imbedded Markov chain. We now wish to deter- 
mine its limiting distribution {a n }^L . The limiting probabilities a n satisfy the equilibrium 
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equations 

«o = aoPo,o + aiPi,o + a 2 P2,o H 

oo 

= 53 a ^:0 (8.2) 

i=0 

&n = Gn-lA) + a nA + O-n+1^2 + " " " 
oo 

= ^a n _i +i /3i, n=l,2, ... (8.3) 

To find the solution of the equilibrium equations it appears that the generating function 
approach does not work here (verify). Instead we adopt the direct approach by trying to 
find solutions of the form 

a n = a n , n = 0,l,2,... (8.4) 

Substitution of this form into equation (8.3) and dividing by the common power a n ~ l yields 



oo 



t=0 



Of course we know that j3i is given by (8.1). Hence we have 

(fXtf 



a = Y,*' T y¥-e-i*f A (t)dt 
£o Jt =° %l 

Jt=o 



The last integral can be recognised as the Laplace-Stieltjes transform of the interarrival 
time. Thus we arrive at the following equation 

a = A{ix- h<j). (8.5) 

We immediately see that a = 1 is a root of equation (8.5), since -4(0) = 1. But this root is 
not useful, because we must be able to normalize the solution of the equilibrium equations. 
It can be shown that (see exercise 50) as long as p < 1 equation (8.5) has a unique root 
a in the range < a < 1, and this is the root which we seek. Note that the remaining 
equilibrium equation (8.2) is also satisfied by (8.4) since the equilibrium equations are 
dependent. We finally have to normalize solution (8.4) yielding 

a n = (1 - a )a n , n = 0,l,2,... (8.6) 

Thus we can conclude that the queue length distribution found just before an arriving 
customer is geometric with parameter a, where a is the unique root of equation (8.5) in 
the interval (0, 1). 
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Example 8.1.1 (M/M/l) 

For exponentially distributed interarrival times we have 

A(s) X 



X + s 

Hence equation (8.5) reduces to 
A 



a 



A + p — pa 
so 

a(X + p — pa) — A = (a — 1)(A — pa) = 0. 

Thus the desired root is a = p and the arrival distribution is given by 

a n = (l- p)p n , n = 0,1,2,... 

Note that this distribution is exactly the same as the equilibrium distribution of the 
M/M/l; see (4.6). This is of course no surprise, because here we have Poisson arrivals. 

Example 8.1.2 (E 2 /M/l) 

Suppose that the interarrival times are Erlang-2 distributed with mean 2/3, so 

Further assume that p = 4 (so p = 3/2 • 1/4 = 3/8 < 1). Then equation (8.5) reduces to 

\7-4aJ 
Thus 

ct(7 - 4a) 2 - 9 = (a - l)(4a - 9)(4a - 1) = 0. 

Hence the desired root is a = 1/4 and 
3 (V 



1U , n = 0,l,2,... 

Example 8.1.3 

Suppose that the interarrival time consist of two exponential phases, the first phase with 
parameter p and the second one with parameter 2p (so it is slightly more complicated than 
Erlang-2 where both phases have the same parameter), where p is also the parameter of 
the exponential service time. The Laplace- Stieltjes transform of the interarrival time is 
given by 

2/i 2 



A(s) 



{p + s)(2p + s) ' 
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For this transform equation (8.5) reduces to 

V _ 2 

{2p — //cr)(3// — fia) (2 — a) (3 — a) 

This leads directly to 

a 3 - 5a 2 + 6a - 2 = (a - l)(a - 2 - V^2)(<r - 2 + y/2) = 0. 
Clearly only the root a = 2 — v2 is acceptable. Therefore we have 

a n = (\/2-i)(2-\/2) n , n = 0,l,2,... 

8.2 Distribution of the sojourn time 

Since the arrival distribution is geometric, it is easy to determine the distribution of the 
sojourn time. In fact, the analysis is similar to the one for for the M/M/l queue (see 
section 4.4). With probability a n an arriving customer finds n customers in the system. 
Then his sojourn time is the sum of n + 1 exponentially distributed service times, each 
with mean \j\i. Hence, 

S(s) = E(e- sS ) 

= £(i-*KpM n+1 



H(l-a 


\ oo 


/i + s 

Mi- 


n=0 



//(l - er) + s ' 

From this we can conclude that the sojourn time S is exponentially distributed with pa- 
rameter //(l — a), i.e., 

P{S <t) = l- e"^ 1 "^, £ > 0. 

Clearly the sojourn time distribution for the G/M/l is of the same form as for the M/M/l, 
the only difference being that p is replaced by a . 

Along the same lines it can be shown that the distributon of the waiting time W is 
given by (cf. (4.14)) 

P{W <t) = l- ae'^ 1 -^, t > 0. 

Note that the probability that a customer does not have to wait is given by 1 — a (and not 
by 1 - p). 
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8.3 Mean sojourn time 

It is tempting to determine the mean sojourn time directly by the mean value approach. 
For an arriving customer we have 

E(S) = E{L a )- + - , (8.7) 

where the random variable L a denotes the number of customers in the system found on 
arrival. According to Little's law it holds that 

E(L) = XE(S). (8.8) 

Unfortunately, we do not have Poisson arrivals, so 

E(L a ) ± E(L). 

Hence the mean value approach does not work here, since we end up with only two equations 
for three unknowns. Additional information is needed in the form of (8.6), yielding 



E{L a ) = Y j na n = Y j n{\ - a)a n = ^~ 



a(l - o)o ,b = -° 

n=0 n=0 

Then it follows from (8.7) and (8.8) that 



E(S) = * +-= - E(L) 



(l-a)fi n (1-a)//' ^'"' (l-o-)yU (1-cr)' 

8.4 Java applet 

There is a JAVA applet avalaible for the evaluation of the G/M/l queue for several distribu- 
tions of the interarrival times. In fact, the applet evaluates the G/E r /1 queue. This queue 
has been analyzed in, e.g., [2]. The link to the applet is http://www.win.tue.nl/cow/Q2. 
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8.5 Exercises 

Exercise 50. 

We want to show that as long as p < 1 equation (8.5) has a unique root a in the range 

< a < 1. Set 

f(a) = A(jji - fia). 

(i) Prove that f(a) is stricly convex on the interval [0, 1], i.e., its derivative is increasing 
on this interval. 

(ii) Show that /(0) > and that /'(l) > 1 provided p < 1. 

(iii) Show that (i) and (ii) imply that the equation a = f(a) has exactly one root in the 
interval (0, 1). 

Exercise 51. 

In a record shop customers arrive according to a hyperexponential arrival process. The 
interarrival time is with probability 1/3 exponentially distributed with a mean of 1 minute 
and with probability 2/3 it is exponentially distributed with a mean of 3 minutes. The 
service times are exponentially distributed with a mean of 1 minute. 

(i) Calculate the distribution of the number of customers found in the record shop by 
an arriving customer. 

(ii) Calculate the mean number of customers in the record shop found on arrival. 

(iii) Determine S(s). 

(iv) Determine the mean time a customer spends in the record shop. 

(v) Calculate the mean number of customers in the record shop (now at an arbitrary 
point in time). 

Exercise 52. 

The distribution of the interarrival time is given by 

F-M = 1 (i " e- s ') + jj (i - «"») . t > 0. 

The service times are exponentially distributed with a mean of 1/6. 

(i) Determine the distribution of the number of customers in the system just before an 
arrival. 

(ii) Determine the distribution of the waiting time. 
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Exercise 53. 

Determine the distribution of the sojourn time in case of exponentially distributed service 

times with mean 1 and hyperexponentially distributed interarrival times with distribution 

function 

F A (t) = 1(1- e-/ 2 ) + \ (l - e-^) , t > 0. 

Exercise 54. 

Consider a queueing system where the interarrival times are exactly 4 minutes. The service 

times are exponentially distributed with a mean of 2 minutes. 

(i) Compute a. 

(ii) Determine the distribution of the sojourn time. 

Exercise 55. 

At a small river cars are brought from the left side to the right side of the river by a ferry. 
On average 15 cars per hour arrive according to a Poisson process. It takes the ferry an 
exponentially distributed time with a mean of 3 minutes to cross the river and return. The 
capacity of the ferry is equal to 2 cars. The ferry only takes off when there are two or more 
cars waiting. 

(i) What is the fraction of time that the ferry is on its way between the two river sides? 

(ii) Determine the distribution of the number of cars that are waiting for the ferry. 

(iii) Determine the mean waiting time of a car. 
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Chapter 9 



Priorities 



In this chapter we analyse queueing models with different types of customers, where one 
or more types of customers have priority over other types. More precisely we consider 
an M/G/l queue with r types of customers. The type i customers arrive according to a 
Poisson stream with rate Aj, i = 1, . . . , r. The service time and residual service of a type 
i customer is denoted by Bi and Ri, respectively. The type 1 customers have the highest 
priority, type 2 customers the second highest priority and so on. We consider two kinds of 
priorities. For the non-preemptive priority rule higher priority customers may not interrupt 
the service time of a lower priority customer, but they have to wait till the service time 
of the low priority customer has been completed. For the preemptive-resume priority rule 
interruptions are allowed and after the interruption the service time of the lower priority 
customer resumes at the point where it was interrupted. In the following two sections we 
show how the mean waiting times can be found for these two kinds of priorities. 

9.1 Non-preemptive priority 

The mean waiting time of a type i customer is denoted by E(Wi) and E(Lj) is the number 
of type i customers waiting in the queue. Further define pi = XiE(Bi). For the highest 
priority customers it holds that 

E(W{) = E(Ll)E(B 1 ) + Y J p J E(R 3 ). 

According to Little's law we have 

E(L\) = A 1 E(W 1 ) 
Combining the two equations yields 

E(W 1 ) = ^= lPjE{R > ) . (9.1) 

1 -pi 
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The determination of the mean waiting time for the lower priority customers is more 
complicated. Consider type % customers with i > 1. The waiting time of a type i customer 
can be divided in a number of portions. The first portion is the amount of work associated 
with the customer in service and all customers with the same or higher priority present 
in the queue upon his arrival. Call this portion Xi. The second portion, say X 2 , is the 
amount of higher priority work arriving during X\. Subsequently the third portion X% is 
the amount of higher priority work arriving during X2, and so on. A realization of the 
waiting time for a type 2 customer is shown in figure 9.1. The increments of the amount 
of work are the service times of the arriving type 1 customers. 

amount of work 




arrivals of type 1 customers 
Figure 9.1: Realization of the waiting time of a type 2 customer 

Hence the mean waiting time is given by 

00 
E(Wi) = E{X X + X 2 + X 3 + •••) = E E(*k)- 

fc=i 

As mentioned above, the first portion of work an arriving type i customer has to wait for 
is the sum of the service times of all customers with the same or higher priority present in 
the queue plus the remaining service time of the customer in service. So 

E{X X ) = JTE{L' 1 j )E{B ] ) + JT Pj E{R 3 ). 
i=i 3=1 

To determine E(X k+ i) note that X k+ i depends on X k . We therefore condition on the 
length of X k . Denote the density of X k by f k (x). Then it follows that 

E(X k+1 ) = / E(X k+l \X k = x)f k {x)dx 

Jx=0 
roo 

= / {XxxE(B 1 ) + ■■■ + \ i ^xE{B i - X ))f k {x)dx 
= (pi + --- + Pi-i)E(X k ). 
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Repeated application of the relation above yields 

E(X k+1 ) = (pi + • • • + Pi . 1 ) k E(X 1 ), k = 0, 1, 2, . . . 
Hence we find for % = 2, . . . , r 

Em = i-( Pl + ... + Pi _ 1 ) i- (pi + . .. + *_,) ■ ^ 9 - 2 ) 

An intuitive argument (which can be made rigorous) to directly obtain the above equation 
is by observing that the waiting time of a type i customer is equal to the first portion of 
work plus all the higher priority work arriving during his waiting time. So 

i-l 

E{Wi) = E(X 1 ) + £ XjEiW^EiBj), 

from which equation (9.2) immediately follows. Substitution of Little's law 

E{Lj) = XiE(Wi) 
into equation (9.2) yields 

{l-{ Pl + ... +Pl ))E{W i ) = 52E(L!)E(B j ) + y EP3E(R3) 

= (l-( Pl + ...+p z . 2 ))E(W l . 1 ). 

By multiplying both sides of this equality with 1 — (pi + . . . + Pj_i) we get the simple 
recursive relation 

(i - X»(! - E^^W) = (i - Eft)(i - EPi)^(Wi-i). 

Repeatedly applying this relation and using (9.1) finally leads to 

E(W t ) = - -. ZUiPjEiRj) i = l,..., r . (9.3) 

1 ; (l-(p 1 + ...+p,))(l-(p 1 + ... + p i _ 1 ))' K ' 

The mean sojourn time E(Si) of a type i customer follows from E(Si) = E(Wi) + E{B i ). i 
yielding 

E M = flvf ( f\ 7 fi + m) > (9-4) 

(1 - (pi + . . . + pi))(l - (pi + . . . + pi_i)) 

for i = 1, • • • , r. 
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9.2 Preemptive-resume priority 

We will show that the results in case the service times may be interrupted easily follow 
from the ones in the previous section. 

Consider a type i customer. For a type i customer there do not exist lower priority 
customers due to the preemption rule. So we henceforth assume that Aj+i = • • • = A r = 0. 

The waiting time of a type i customer can again be divided into portions X\,X 2 , .... 
Now Xi is equal to the total amount of work in the system upon arrival, since we assumed 
that there are no lower priority customers. Observe that the total amount of work in 
the system does not depend on the order in which the customers are served. Hence, at 
each point in time, it is exactly the same as in the system where the customers are served 
according to the non-preemptive priority rule. So Xi,X 2 , . . ., and thus also Wi have the 
same distribution as in the system with non-preemptive priorities and, of course, with 
Aj + i = • • • = A r = 0. From (9.3) we then obtain 

E<W t )= gkftgW i = 1 ... r 

For the mean sojourn time we have to add the service time plus all the interruptions of 
higher priority customers during the service time. The mean of such a generalized service 
time can be found along the same lines as (9.2), yielding 

E(Bj) 

1 - (Pi H h Pi-i) ' 

So the mean sojourn time of a type % customer is given by 

E(s) = YUPiEjRi) | gW (Q5) 

1 l) (l-(p 1 + ...+A))(l-(p 1 + ... + p i _ 1 )) + l-(p 1 + ---+p,_ 1 )' l ' } 

for i — 1, • • • , r. 

9.3 Shortest processing time first 

In production systems one often processes jobs according to the shortest processing time 
first rule (SPTF). The mean production lead time in a single machine system operating 
according to the SPTF rule can be found using the results in section 9.1. 

Consider an M/G/l queue with arrival rate A and service times B with density /b(x). 
Assume that p = XE(B) < 1. The server works according to the SPTF rule. That is, after 
a service completion, the next customer to be served is the one with the shortest service 
time. 

Define type x customers as the ones with a service time between x and x + dx. The 
mean waiting time of a type x customer is denoted by E(W(x)) and p(x)dx is the fraction 
of time the server helps type x customers, so 

p(x)dx = (XfB(x)dx)x = XxfB(x)dx. (9.6) 
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From (9.3) and by observing that the numerator in (9.3) corresponds to the mean amount 
of work at the server, which in the present situation is simply given by pE(R), we obtain 

E(W(x)) . pE ^ 



(i - I y % x p{v)dyf 

pE(R) 

(i-Xf y % x yfB(y)d y y 



Hence the mean overall waiting time is given by 

poo 

E(W) = / E(W(x))f B (x)dx 

Jx=0 



nF /m r fB(x)dx 

= P E \R) / T. x ry=x t , \j \ 2 ■ 9 - 7 

Jx=o (1 - A/JL,, yf B (y)dy) 2 

In table 9.1 we compare the SPTF rule with the usual first come first served (FCFS) rule 
for an M/M/l system with mean service time 1. The mean waiting time for the SPTF 
rule is given by 



roo 

E(W) = p / 

Jx=0 1 — 0(1 



e x dx 



p(l — e x — xe x )) 2 

and for the FCFS rule it satisfies 
P 



E(W) 



1-p 



The results in table 9.1 show that considerable reductions in the mean waiting time are 
possible. 





E 


\W) 


p 


FCFS 


SPTF 


0.5 


1 


0.713 


0.8 


4 


1.883 


0.9 


9 


3.198 


0.95 


19 


5.265 



Table 9.1: The mean waiting time for FCFS and SPTF in an M/M/l with E(B) 



9.4 A conservation law 

In this section we consider a single-server queue with r types of customers. The type % 
customers arrive according to a general arrival stream with rate Aj, i = 1, . . . , r. The mean 
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service time and mean residual service time of a type i customer is denoted by E(B{) and 
E(Ri), respectively. Define pi = XiE(Bi). We assume that 

r 

£A<£(B 4 )<1, 

i=l 

so that the server can handle the amount of work offered per unit of time. Customers enter 
service in an order independent of their service times and they may not be interrupted 
during their service. So, for example, the customers may be served according to FCFS, 
random or a non-preemptive priority rule. Below we derive a conservation law for the 
mean waiting times of the r types of customers, which expresses that a weighted sum of 
these mean waiting times is independent of the service discipline. This implies that an 
improvement in the mean waiting of one customer type owing to a service discipline will 
always degrade the mean waiting time of another customer type. 

Let E{V{P)) and E(Lf(P)) denote the mean amount of work in the system and the 
mean number of type i customers waiting in the queue, respectively, for discipline P. The 
mean amount of work in the system is given by 

E(V(P)) =±E(Lj(P))E(B. l ) + ± Pi E(R l ). (9.8) 

i=l »=1 

The first sum at the right-hand side is the mean amount of work in the queue, and the 
second one is the mean amount of work at the server. Clearly the latter does not depend 
on the discipline P. 

The crucial observation is that the amount of work in the system does not depend 
on the order in which the customers are served. The amount of work decreases with one 
unit per unit of time independent of the customer being served and when a new customer 
arrives the amount of work is increased by the service time of the new customer. Hence, 
the amount of work does not depend on P. Thus from equation (9.8) and Little's law 

E{L\) = \EiWiiP)), 

we obtain the following conservation law for the mean waiting times, 

r 

y^ pjE(Wj(P)) = constant with respect to service discipline P. 
%=\ 

Below we present two examples where this law is used. 

Example 9.4.1 (M/G/l with FCFS and SPTF) 

The SPTF rule selects customers in a way that is dependent of their service times. Nev- 
ertheless, the law above also applies to this rule. The reason is that the SPTF rule can be 
translated into a non-preemptive rule as explained in section 9.3. Below we check whether 
for the M/G/l the weighted sum of the mean waiting times for SPTF is indeed the same 
as for FCFS. 
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In case the customers are served in order of arrival it holds that (see (7.14)) 

1 - p 

When the server works according to the SPTF rule we have (see (9.6) and (9.7)) 

pE (R) Xx fB(x)dx 



E(W(x))p(x)dx = = , 

*=o Jx=o (1- A/jLo yjBWdyy 

P E{R) 



x=0 



i-^I y y : x yfB(y)dy 
P 2 E(R) 

1-p ' 

which indeed is the same as for the FCFS rule. 

Example 9.4.2 (M/G/l with non-preemptive priority) 

Consider an M/G/l queue with two types of customers. The type 1 customers have 
non-preemptive priority over the type 2 customers. The mean waiting time for the type 1 
customers is given by (9.1). We now derive the mean waiting time for the type 2 customers 
by using the conservation law. According to this law it holds that 

pxEiWt) + p 2 E(W 2 ) = C, (9.9) 

where C is some constant independent of the service discipline. To determine C consider 
the FCFS discipline. For FCFS it follows from (7.14) that 



1 - pi - p 2 
Hence, 

1 - p x - p 2 

By substituting (9.1) and (9.10) into equation (9.9) we retrieve formula (9.3) for the mean 
waiting time of the type 2 customers under the non-preemptive priority rule. 
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9.5 Exercises 

Exercise 56. 

Customers arrive to a single-server queue according to a Poisson process with a rate of 10 
customers per hour. Half of the customers has an exponential service time with a mean of 
3 minutes. The other half has an exponential service time with a mean of 6 minutes. The 
first half are called type 1 customers, the second half type 2 customers. Determine in each 
of the following three cases the mean waiting time: 

(i) Service in order of arrival (FCFS). 
(ii) Non-preemptive priority for type 1 customers, 
(ii) Non-preemptive priority for type 2 customers. 

Exercise 57. 

Consider an M/D/l queue with three types of customers. The customers arrive according 

to a Poisson process. Per hour arrive on average 1 type 1 customer, 2 type 2 customers 

and also 2 type 3 customers. Type 1 customers have preemptive resume priority over 

type 2 and 3 customers and type 2 customers have preemptive resume priority over type 3 

customers. The service time of each customer is 10 minutes. 

Calculate the mean sojourn time for each of the three types of customers. 

Exercise 58. 

Consider a machine where jobs arrive according to a Poisson stream with a rate of 4 jobs 
per hour. Half of the jobs have a processing time of exactly 10 minutes, a quarter of 
the jobs have a processing time of exactly 15 minutes and the remaining quarter have a 
processing time of 20 minutes. The jobs with a processing time of 10 minutes are called 
type 1 jobs, the ones with a processing time of 15 minutes type 2 jobs and the rest type 3 
jobs. The jobs are processed in order of arrival. 

(i) Determine the mean sojourn time (waiting time plus processing time) of a type 1, 2 
and 3 job and also of an arbitrary job. 

One decides to process smaller jobs with priority. So type 1 orders have highest priority, 
type 2 orders second highest priority and type 3 orders lowest priority. 
Answer question (i) for the following two cases: 

(ii) Jobs processed at the machine may not be interrupted. 

(iii) Type 1 and type 2 jobs may interrupt the processing of a type 3 job. Type 1 jobs 
may not interrupt the processing of a type 2 job. 

Exercise 59. 

A machine produces a specific part type. Orders for the production of these parts arrive 
according to a Poisson stream with a rate of 10 orders per hour. The number of parts 
that has to be produced for an order is equal to n with probability 0.5 n , n = 1, 2, . . . The 
production of one part takes exactly 2 minutes. Orders are processed in order of arrival. 
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(i) Denote by B the production time in minutes of an arbitrary order. Show that 

E(B)=4, a 2 (B) = 8. 

(ii) Determine the mean production lead time (waiting time plus production time) of an 
arbitrary order. 

The management decides to give orders for the production of 1 part priority over the other 
orders. But the production of an order may not be interrupted. 

(iii) Determine the mean production lead of an order for the production of 1 part, and of 
an order for the production of at least 2 parts. 

(iv) Determine the mean production lead time of an arbitrary order. 

Exercise 60. 

Consider a machine for mounting electronic components on printed circuit boards. Per 
hour arrive according to a Poisson process on average 30 printed circuit boards. The 
time required to mount all components on a printed circuit board is uniformly distributed 
between 1 and 2 minutes. 

(i) Calculate the mean sojourn time of an arbitrary printed circuit board. 

One decides to give printed circuit boards with a mounting time between 1 and x (1 < 
x < 2) minutes non-preemptive priority over printed circuit boards the mounting time of 
which is greater than x minutes. 

(ii) Determine the mean sojourn time of an arbitrary printed circuit board as a function 
of x. 

(iii) Determine the value of x for which the mean sojourn time of an arbitrary printed 
circuit board is minimized, and calculate for this specific x the relative improvement 
with respect to the mean sojourn time calculated in (i). 

Exercise 61. 

A machine produces 2 types of products, type A and type B. Production orders (for one 
product) arrive according to a Poisson process. For the production of a type A product 
arrive on average 105 orders per day (8 hours), and for a type B product 135 orders per 
day. The processing times are exponentially distributed. The mean processing time for a 
type A order is 1 minute and for a type B order it is 2 minutes. Orders are processed in 
order of arrival. 

(i) Calculate the mean production lead time for a type A product and for a type B 
product. 

(ii) Determine the Laplace-Stieltjes transform of the waiting time. 
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(iii) Determine the Laplace-Stieltjes transform of the production lead time of a type A 
product and determine the distribution of the production lead time of a type A 
product. 

For the production lead time of a type A product one uses a norm of 10 minutes. Each 
time the production lead time of a type A product exceeds this norm, a cost of 100 dollar 
is charged. 

(iv) Calculate the average cost per day. 

To reduce the cost one decides to give type A products preemptive resume priority over 
type B products. 

(v) Determine the mean production lead time for a type A product and for a type B 
product. 

(vi) Calculate the average cost per day. 

Exercise 62. 

A machine mounts electronic components on three different types of printed circuit boards, 
type A, B and C boards say. Per hour arrive on average 60 type A boards, 18 type B 
boards and 48 type C boards. The arrival streams are Poisson. The mounting times are 
exactly 20 seconds for type A, 40 seconds for type B and 30 seconds for type C. The 
boards are processed in order of arrival. 

(i) Calculate for each type of printed circuit board the mean waiting time and also 
calculate the mean overall waiting time. 

Now suppose that the printed circuit boards are processed according to the SPTF rule. 

(ii) Calculate for each type of printed circuit board the mean waiting time and also 
calculate the mean overall waiting time. 
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Chapter 10 

Variations of the M/G/l model 



In this chapter we treat some variations of the M/G/l model and we demonstrate that the 
mean value technique is a powerful technique to evaluate mean performance characteristics 
in these models. 



10.1 Machine with setup times 

Consider a single machine where jobs are being processed in order of arrival and suppose 
that it is expensive to keep the machine in operation while there are no jobs. Therefore the 
machine is turned off as soon as the system is empty. When a new job arrives the machine 
is turned on again, but it takes some setup time till the machine is ready for processing. 
So turning off the machine leads to longer production leadtimes. But how much longer? 
This will be investigated for some simple models in the following subsections. 



10.1.1 Exponential processing and setup times 

Suppose that the jobs arrive according to a Poisson stream with rate A and that the 
processing times are exponentially distributed with mean 1/ p. For stability we have to 
require that p = X/p < 1. The setup time of the machine is also exponentially distributed 
with mean 1/0. We now wish to determine the mean production lead time E(S) and the 
mean number of jobs in the system. These means can be determined by using the mean 
value technique. 

To derive an equation for the mean production lead time, i.e. the arrival relation, we 
evaluate what is seen by an arriving job. We know that the mean number of jobs in the 
system found by an arriving job is equal to E(L) and each of them (also the one being 
processed) has an exponential (residual) processing time with mean 1/ p. With probability 
1 — p the machine is not in operation on arrival, in which case the job also has to wait 
for the (residual) setup phase with mean 1/0. Further the job has to wait for its own 
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processing time. Hence 

E(S) = (l-p) 1 - + E(L)- + - 

U p p 

and together with Little's law, 

E(L) = XE(S) 
we immediately find 



E(S) 



l- p e 



So the mean production lead time is equal to the one in the system where the machine is 
always on, plus an extra delay caused by turning off the machine when there is no work. 
In fact, it can be shown that the extra delay is exponentially distributed with mean 1/9. 

10.1.2 General processing and setup times 

We now consider the model with generally distributed processing times and generally 
distributed setup times. The arrivals are still Poisson with rate A. The first and second 
moment of the processing time are denoted by E(B) and E(B 2 ) respectively, E(T) and 
E(T 2 ) are the first and second moment of the setup time. For stability we require that 
p = XE(B) < 1. Below we demonstrate that also in this more general setting the mean 
value technique can be used to find the mean production lead time. We first determine the 
mean waiting time. Then the mean production lead time is found afterwards by adding 
the mean processing time. 

The mean waiting time E(W) of a job satisfies 

E(W) = E(L q )E(B) + pE(R B ) 

+ P(Machine is off on arrival)i?(T) 

+ P(Machine is in setup phase on amval) E (Rt) , (10.1) 

where E(Rb) and E(R T ) denote the mean residual processing and residual setup time, so 
(see (7.15)) 

F(R ) - E{B2) F(R ) - EiT2) 

e{Rb) ~ 2E\B) ' ( t) ~ mf) ' 

To find the probability that on arrival the machine is off (i.e. not working and not in the 
setup phase), note that by PASTA, this probability is equal to the fraction of time that 
the machine is off. Since a period in which the machine is not processing jobs consists of 
an interarrival time followed by a setup time, we have 

P (Machine is off on arrival) = (1 — p) — ; — -. 

v ' v HJ l/\ + E(T) 
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Similarly we find 

E(T) 



P(Machine is in setup phase on arrival) = (1 — p)- 

Substituting these relations into (10.1) and using Little's law stating that 

E(L q ) = \E(W) 
we finally obtain that 

EM) = eEM + Uh BIT) + E{T) E(R T ) 

Note that the first term at the right-hand side is equal to the mean waiting time in the 
M/G/l without setup times, the other terms express the extra delay due to the setup 
times. Finally, the mean production lead time E(S) follows by simply adding the mean 
processing time E(B) to E(W). 

10.1.3 Threshold setup policy 

A natural extension to the setup policy in the previous section is the one in which the 
machine is switched on when the number of jobs in the system reaches some threshold 
value, TV say. This situation can be analyzed along the same lines as in the previous 
section. The mean waiting time now satisfies 

E(W) = E(L q )E(B) + pE(R B ) 

v^ (N-i \ 

+ 2j -P(Arriving job is number i) ( — \- E[T) J 

+ P(Machine is in setup phase on arrival)P(i?T), (10-2) 

The probability that an arriving job is the z-th one in a cycle can be determined as follows. 
A typical production cycle is displayed in figure 10.1. 

i^ Non-processing period ^j^ Processing period ^j 

A, A 9 ... A N Setup B, B 9 

I — II h-^H 1 — -^-^ 1 1 

Figure 10.1: Production cycle in case the machine is switched on when there are N jobs 

The probability that a job arrives in a non-processing period is equal to 1 — p. Such a 
period now consists of iV interarrival times followed by a setup time. Hence, the probability 
that a job is the i-th one, given that the job arrives in a non-processing period, is equal to 
1/A divided by the mean length of a non-processing period. So 

P(Arriving job is number z) = (1 - p) — — i = l,...,N, 
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and similarly, 

P (Machine is in setup phase on arrival) = (1 — p) — — -. 

Substituting these relations into (10.2) we obtain, together with Little's law, that 

f(w\- pE{Rb) \ N,X ( N - 1 + f(t)\a e{t) F(R ) 

E{W) - T^T + N/X + E(T) \-2X~ + E{T) ) + N/X + E(T) E{Rt) - 

10.2 Unreliable machine 

In this section we consider an unreliable machine processing jobs. The machine can break 
down at any time it is operational, even though it is not processing jobs. What is the impact 
of these breakdowns on the production leadtimes? To obtain some insight in the effects of 
breakdowns we formulate and study some simple models in the following subsections. 

10.2.1 Exponential processing and down times 

Suppose that jobs arrive according to a Poisson stream with rate A. The processing times 
are exponentially distributed with mean 1/ ' \x. The machine is successively up and down. 
The up and down times of the machine are also exponentially distributed with means 1/rj 
and 1/6, respectively. 

We begin with formulating the condition under which the machine can handle the 
amount of work offered per unit time. Let pjj and po denote the fraction of time the 
machine is up and down, respectively. So 

1/rj 1 1 



PU = T~, T-TfT, = ^—, 77; » PD = 1 - Pu 



I/77 + I/0 l + rj/6' ru ru 1 + 6/r]' 

Then we have to require that 

- < Pu- (10.3) 

P 

We now proceed to derive an equation for the mean production lead time. An arriving 
job finds on average E(L) jobs in the system and each of them has an exponential processing 
time with mean 1/p. So in case of a perfect machine his mean sojourn time is equal to 
(E(L) + l)/p. However, it is not perfect. Breakdowns occur according to a Poisson 
process with rate 77. So the mean number of breakdowns experienced by our job is equal to 
r](E(L) + l)/p, and the mean duration of each breakdown is 1/6. Finally, with probability 
pD the machine is already down on arrival, in which case our job has an extra mean delay 
of 1/6. Summarizing we have 

E(S) = (E(L) + l)- + r,(E(L) + l)-.l + p D - 

p p V V 
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Then, with Little's law stating that 

E(L) = XE(S), 

we immediately obtain 

E(s) = ^KvPu)+ Pd/0 
1 ' l-\/(jipu) 

In table 10.1 we investigate the impact of the variability in the availability of the 
machine on the mean production leadtime. The mean production time is 1 hour (// = 1). 
The average number of jobs that arrive in a week (40 hours) is 32, so A = 0.8 jobs per 
hour. The fraction of time the machine is available is kept constant at 90%, so pu = 0.9 
(and thus r]/9 = 1/9). The rate r\ at which the machine breaks down is varied from (on 
the average) every 10 minutes till once a week. In the former case the mean down time is 
1.1 minute, in the latter case it is more dramatic, namely nearly half a day (4.4 hours). 
The results show that the variation in the availability of the machine is essential to the 
behavior of the system. Note that as rj and 9 both tend to infinity such that rj/6 = 1/9, 
then E(S) tends to 10, which is the mean sojourn time in an M/M/l with arrival rate 0.8 
and service rate 0.9. 



V 


9 


E(S) 


6 


54 


10.02 


3 


27 


10.03 


1 


9 


10.1 


0.125 


1.125 


10.8 


0.0625 


0.5625 


11.6 


0.025 


0.225 


14 



Table 10.1: The mean production leadtime in hours as a function of the break-down rate 
rj per hour for fixed availability of 90% 



10.2.2 General processing and down times 

We now consider the model with general processing times and general down times. The 
arrivals are still Poisson with rate A. The first and second moment of the processing 
time are denoted by E(B) and E(B 2 ) respectively. The time between two breakdowns is 
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exponentially distributed with mean 1/r] and E(D) and E(D 2 ) are the first and second 
moment of the down time. For stability we require that (cf. (10.3)) 

XE(B) < ^tttt. 

v ; l + r)E(D) 

Below we first determine the mean waiting time by using the mean value technique. The 
mean production leadtime is determined afterwards. 

We start with introducing the generalized processing time, which is defined as the pro- 
cessing time plus the down times occurring in that processing time. Denote the generalized 
processing time by G. Then we have 

N(B) 

8=1 

where N(B) is the number of break-downs during the processing time B and Dj is the z-th 
down time. For the mean of G we get by conditioning on B and N(B) that 

.(v x Y 

Cj\^r\n = x, iv \n ) = ri)t • 
\r]x) n 



poo ^ \nr) 

E(G) = / T E(G\B = x } N(B)=n)e- llxK - L ^-f B (x)dx 
Jx=o ^ n! 

r°° °° <nr) n 

= / y £( x + nE(D))e- rpi ^-fB(x)dx 



'X; 



n=0 

■DC' 



J2(x + x V E(D))e- r ' x ^f B (x)dx 



x=o n=0 n 

= E(B) + E(B)rjE(D), 

and similarly, 

E(G 2 ) = E(B 2 )(1 + T)E(D)) 2 + E(B)riE(D 2 ). 

A typical production cycle is shown in figure 10.2. The non-processing period consists of 
cycles which start with an up time. When the machine is up two things can happen: a job 
arrives, in which case the machine starts to work, or the machine goes down and has to 
be repaired. Hence, an up time is exponentially distributed with mean 1/(A + rj) and it is 
followed by a processing period (i.e. a job arrives) with probability \/(\ + r)) or otherwise, 
it is followed by a down time. During a processing period the machine works on jobs with 
generalized processing times (until all jobs are cleared). For the mean waiting time it holds 
that 

E(W) = E{L*)E{G) + p G E(R G ) 

+P[Arrival in a down time in a non-processing period] E(Rd), (10.4) 

where p G is the fraction of time the machine works on generalized jobs and E(Rq) and 
E(Rd) denote the mean residual generalized processing time and the mean residual down 
time, so 

p G = XE(G), E(R G ) = ^)-,E(R D ) Lil) " ] 



2E(G) ' v ' 2E{D) ' 
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Ug Non-processing period ^l^ Processing period ^i 

up down up down Gi G9 

I 1 1 1 1 1 1 1 1 



first arrival 
Figure 10.2: Production cycle of an unreliable machine 

The probability that a job arriving in a non-processing period finds that the machine is 
down, is given by 

E(D) V /{X + V ) EiD)^ 



1/{X + V ) + E{D) V /{X + V ) l + E(D)ri' 

Hence, since a job arrives with probability 1 — pa in a non-processing period, we have 

FY Tl^n 
P(Arrival in a down time in a non-processing period) = (1 — pa) . 

Substitution of this relation into (10.4) and using Little's law yields 

and the mean production lead time finally follows from 

E(S) = E(W) + E{G). 

10.3 M/G/l queue with an exceptional first customer 
in a busy period 

A servers' life in an M/G/l queue is an alternating sequence of periods during which 
no work is done, the so-called idle periods, and periods during which the server helps 
customers, the so-called busy periods. In some applications the service time of the first 
customer in a busy period is different from the service times of the other customers served 
in the busy period. For example, consider the setup problem in section 10.1 again. This 
problem may, alternatively, be formulated as an M/G/l queue in which the first job in 
a busy period has an exceptional service time, namely the setup time plus his actual 
processing time. The problem in section 10.2 can also be described in this way. Here we 
can take the generalized processing times as the service times. However, on arrival of the 
first job in a busy period the machine can be down. Then the (generalized) servicing of 
that job cannot immediately start, but one has to wait till the machine is repaired. This 
repair time can be included in the service time of the first job. In this case, however, the 
determination of the distribution, or the moments of the service time of the first job is 
more difficult then in the setup problem. Below we show how, in the general setting of an 
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M/G/l with an exceptional first customer, the mean sojourn time can be found by using 
the mean- value approach. 

Let Bf and Rf denote the service time and residual service time, respectively, of the 
first customers in a busy period. For the mean waiting time we have 

E(W) = E{L q )E(B) + pfE(R f ) + pE(R), (10.5) 

where pf is the fraction of time the server works on first customers, and p is the fraction of 
time the server works on other (ordinary) customers. Together with Little's law we then 
obtain from (10.5) that 

and for the mean sojourn time we get 

E(S) = E(W) + (1 - Pf - p)E(B f ) + (p f + p)E(B). 

It remains to determine pf and p. The probability that an arriving customer is the first 
one in a busy cycle is given by 1 — pf — p. Hence, the number of first customers arriving 
per unit of time is equal to A(l — pf — p). Thus pf and p satisfy 

Pf = X(l-p f -p)E(Bf) 1 
p = X(p f + p)E(B), 

from which it follows that 

XE(Bf)(l-XE(B)) _ XE(B f )XE(B) 

Pf ~ 1 + XE(B f ) - XE(B) ' P ~ l + XE(B f )-XE(B)' 

10.4 M/G/l queue with group arrivals 

In this section we consider the M/G/l queue where customers do not arrive one by one, 
but in groups. These groups arrive according to a Poisson process with rate A. The group 
size is denoted by the random variable G with probability distribution 

g k = P(G = k), fc = 0,l,2,... 

Note that we also admit zero-size groups to arrive. Our interest lies in the mean waiting 
time of a customer, for which we can write down the following equation. 

oo 

EiW) = E(L q )E(B) + pE(R) + ]£ r k (k - l)E(B) , (10.6) 

fc=i 

where p is the server utilization, so 
p = XE(G)E(B), 
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and rfc is the probability that our customer is the fcth customer served in his group. The 
first two terms at the right-hand side of (10.6) correspond to the mean waiting time of 
the whole group. The last one indicates the mean waiting time due to the servicing of 
members in his own group. 

To find ffc we first determine the probability h n that our customer is a member of a 
group of size n (cf. section 7.7). Since it is more likely that our customer belongs to a 
large group than to a small one, it follows that h n is proportional to the group size n as 
well as the frequency of such groups. Thus we can write 

h n = Cng n , 

where C is a constant to normalize this distribution. So 

oo 

C- 1 = J2ng n = E(G). 

n=l 

Hence 

, _ ng n _ 

n ~ E(G)> "- i > / >--- 

Given that our customer is a member of a group of size n, he will be with probability \/n 
the fcth customer in his group going into service (of course, n> k). So we obtain 

OO 1 1 oo 

r k = Z^ h n ■ — = JJJTqs 2^ 9n, 

n=k \ I n=k 

and for the last term in (10.6) it immediately follows that 

OO -I oo oo 

Y.r k {k-l)E{B) = -^Y,Y,9n(h-l)E{B) 

fe=l ^V^) k=ln=k 

-i oo n 

M G ) n=lk=l 
1 °° 1 

J2^n(n-l)g n E(B) 



E(G) £i 2 
E{G 2 ) - E(G) 



E{B). (10.7) 



2E(G) 
From (10.6) and (10.7) and Little's law stating that 

E(L q ) = XE{G)E(W), 
we finally obtain 

man pE(H) j. (g(G 2 ) - E(G))E(B) 
E{W) = T^ + 2B(0(l- rt ' 

The first term at the right-hand side is equal to the mean waiting time in the system where 
customers arrive one by one according to a Poisson process with rate \E(G). Clearly, the 
second term indicates the extra mean delay due to clustering of arrivals. 
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Example 10.4.1 (Uniform group sizes) 

In case the group size is uniformly distributed over 1,2, ... ,n, so 

9k = - , k=l,...,n, 

n 



we find 



E(G) = ± k - = n -±l, E^-G) = ± k2 - k -^-^^) 



k =i n z S n 



Hence 

1 ] 1-P 3(1 -p) 
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10.5 Exercises 

Exercise 63. 

Consider an M/G/l queue where at the end of each busy period the server leaves for a 

fixed period of T units of time. Determine the mean sojourn time E(S). 

Exercise 64. 

In a factory a paternoster elevator (see figure 10.3) is used to transport rectangular bins 
from the ceiling to the floor. This is a chain-elevator with product carriers (each carrier 
can hold at most one bin) at a certain fixed distance. The number of carriers between the 
ceiling and the floor is 4. Every time a carrier reaches the ceiling (and, simultaneously, 
another one reaches the floor) the elevator stops 2 seconds. This is exactly the time a bin 
needs to enter or leave the carrier. Subsequently the elevator starts to move again till the 
next carrier reaches the ceiling. This takes exactly 3 seconds. Then it waits for 2 seconds, 
moves again, and so on. The process of waiting and moving goes on continuously. Note 
that a bin, finding on arrival an empty carrier positioned in front of it, cannot enter the 
carrier, because the residual time till the carrier start to move again is less than 2 seconds. 
The bin has to wait till the next carrier arrives. 



m 



m 



Figure 10.3: Paternoster elevator 

Bins arrive according to a Poisson process with a rate of 10 bins per minute, 
(i) Determine the mean waiting time of a bin in front of the elevator. 

(ii) Determine the mean time that elapses from the arrival of a bin till the moment the 
bin leaves the elevator downstairs. 

Exercise 65. 

Consider a machine which is turned off when there is no work. It is turned on and restarts 
work when enough orders, say N, arrived to the machine. The setup times are negligible. 
The processing times are exponentially distributed with mean 1/ '// and the average number 
of orders arriving per unit time is A(< fi). 
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(i) Determine the mean number of orders in the system. 

(ii) Determine the mean production lead time. 

Suppose that A = 20 orders per hour, 1/fi = 2 minutes and that the setup cost is 54 
dollar. In operation the machine costs 12 dollar per minute. The waiting cost is 1 dollar 
per minute per order. 

(iii) Determine, for given threshold N, the average cost per unit time. 

(iv) Compute the threshold N minimizing the average cost per unit time. 

Exercise 66. 

At a small river pedestrians are brought from the left side to the right side of the river 

by a ferry. On average 40 pedestrians per hour arrive according to a Poisson process. It 

takes the ferry exactly 2 minutes to cross the river and return. The capacity of the ferry 

is sufficiently big, so it is always possible to take all waiting pedestrians to the other side 

of the river. The ferry travels continuously back and forth, also when there are no waiting 

pedestrians. 

(i) Determine the mean waiting time and the mean number of pedestrians waiting for 
the ferry. 

The ferry now only takes off when there are one or more pedestrians. In case there are no 
pedestrians the ferry waits for the first one to arrive, af ter which it immediately takes off 
(the ferry never waits at the other side of the river). 

(ii) Determine the probability that an arriving pedestrian finds the ferry waiting to take 
off. 

(iii) Determine the mean waiting time and the mean number of pedestrians waiting for 
the ferry. 

Exercise 67. 

A machine produces products in two phases. The first phase is standard and the same for 
all products. The second phase is customer specific (the finishing touch). The first (resp. 
second) phase takes an exponential time with a mean of 10 (resp. 2) minutes. Orders for 
the production of one product arrive according to a Poisson stream with a rate of 3 orders 
per hour. Orders are processed in order of arrival. 

(i) Determine the mean production lead time (waiting time plus production time) of an 
order. 

The machine is switched off when the system is empty and it is switched on again as soon 
as the first order arrives. A fixed cost of 20$ is incurred each time the machine is switched 
on (the time needed to switch the machine on or off is negligible). 
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(ii) Determine the average switch-on cost per hour. 

To reduce the production lead time one decides to start already with the production of 
phase 1 when the system is empty. If upon completion of phase 1 no order has arrived 
yet, the production stops and the machine is switched off. When the first order arrives the 
machine is switched on again and can directly start with phase 2. 

(iii) Determine the reduction in the mean production lead time. 

(iv) Determine the average switch-on cost per hour. 

Exercise 68. 

In a bank there is a special office for insurances. Here arrive according to a Poisson process 
on average 8 customers per hour. The service times are exponentially distributed with a 
mean of 5 minutes. As soon as all customers are served and the office is empty again, the 
clerck also leaves to get some coffee. He returns after an exponential time with a mean of 
5 minutes and then starts servicing the waiting customers (if there are any) or patiently 
waits for the first customer to arrive. 

(i) What is the mean sojourn time and the mean number of customers in the office? 

Now suppose that when the clerck finds an empty office upon his return, he immediately 
leaves again to get another cup of coffee. This takes again an exponential time with a mean 
of 5 minutes. The clerck repeats to leave until he finds a waiting customer upon return. 

(ii) How many cups of coffee the clerck drinks on average before he starts servicing again? 

(iii) Determine the mean sojourn time and the mean number of customers in the office. 

Exercise 69. 

Consider a machine processing orders. These orders arrive according to a Poisson process 
with a rate of 1 order per hour. When there are no orders, the machine is switched off. As 
soon as a new order arrives, the machine is switched on again, which takes exactly T hours. 
The processing time of an order is exponentially distributed with a mean of 30 minutes. 

(i) Determine as a function of T: 

(a) The mean numbers of orders processed in a production cycle. 

(b) The mean duration of a production cycle. 

(c) The mean production lead time of an order. 

Suppose that it costs 17 dollar each time the machine is switched on again and that the 
waiting cost per hour per order are 1 dollar. 

(ii) Show that the average cost per hour are minimal for T = 3 hour. 
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Exercise 70. 

Consider a queueing system where on average 3 groups of customers arrive per hour. The 
mean group size is 10 customers. The service time is exactly 1 minute for each customer. 
Determine the mean sojourn time of the first customer in a group, the last customer in a 
group and an arbitrary one in the following two cases: 

(i) the group size is Poisson distributed; 

(ii) the groups size is geometrically distributed. 

Exercise 71. 

Customers arrive in groups at a server. The groups consist of 1 or 3 customers, both with 
equal probability, and they arrive according to a Poisson stream with a rate of 2 groups 
per hour. Customers are served one by one, and they require an exponential service time 
with a mean of 10 minutes. 

(i) Determine the mean sojourn time of an arbitrary customer. 

(ii) Determine the mean number of customers in the system. 

Exercise 72. 

Passengers are brought with small vans from the airport to hotels nearby. At one of 

those hotels on average 6 vans per hour arrive according to a Poisson process. With 

probability 1/4 a van brings 2 guests for the hotel, with probability 1/4 only one guest 

and with probability 1/2 no guests at all. At the reception of the hotel there is always one 

receptionist present. It takes an exponential time with a mean of 5 minutes to check in a 

guest. 

(i) Determine the distribution of the number of guests at the reception. 

(ii) Determine the mean waiting time of an arbitrary guest at the reception. 

Exercise 73. 

Customers arrive at a server according to a Poisson stream with a rate of 6 customers per 
hour. As soon as there are no customers the server leaves to do something else. He returns 
when there are 2 customers waiting for service again. It takes exactly 5 minutes for him to 
return. The service time of a customer is uniformly distributed between 5 and 10 minutes. 

(i) Determine the mean duration of a busy period, i.e., a period during which the server 
is servicing customers without interruptions. 

(ii) Determine the mean number of customers served in a busy period. 

(iii) Determine the mean sojourn time of a customer. 
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Chapter 11 
Insensitive systems 



In this chapter we study some queueing systems for which the queue length distribution is 
insensitive to the distribution of the service time, but only depends on its mean. 

11.1 M/G/oo queue 

In this model customers arrive according to a Poisson process with rate A. Their service 
times are independent and identically distributed with some general distribution function. 
The number of servers is infinite. So there is always a server available for each arriving 
customer. Hence, the waiting time of each customer is zero and the sojourn time is equal 
to the service time. Thus by Little's law we immediately obtain that 

E(L) = p, 

where p = XE(B) denotes the mean amount of work that arrives per unit time. In the 
remainder of this section we want to also determine the distribution of L, i.e., the proba- 
bilities p n that there are n customers in the system. 

Example 11.1.1 (M/M/oo) 

In this model the service times are exponentially distributed with mean 1/p. 

XX X 





a 2a nu, 

Figure 11.1: Flow diagram for the M/M/oo model 

From figure 11.1 we obtain by equating the flow from state n — 1 to n and the flow 
from n to n — 1 that 

p n -iX = p n np. 

Ill 



Thus 

X p p 2 p n 

Pn = Pn-1 = -Pn-1 = ~, 7TPn-2 = • • • = — p . 

np n n[n — 1) n\ 

Since the probabilities p n have to add up to one, it follows that 

oo n 
-1 _ V^ P _ „p 



po" 1 = £ 



e 



„=o n! 



Summarizing, we have found that 

Pn = -,e~ p . (11.1) 

n! 

Thus the number of customers in the system has a Poisson distribution with mean p. 

Example 11.1.2 (M/D/oo) 

Let b denote the constant service time. The probability p n (t) that there are exactly n 
customers in the system at time t is equal to the probability that between time t — b and t 
exactly n customers arrived. Since the number of customers arriving in a time interval of 
length b is Poisson distributed with mean Xb, we immediately obtain 

Mt) - <Te- = ^, 
n\ n\ 

which is valid for all t > b, and thus also for the limiting distribution. 

Example 11.1.3 (Discrete service time distribution) 

Suppose that the service time distribution is discrete, i.e., there are nonnegative numbers 

bi and probabilities <& such that 

P(B = b l ) = q l , i = l,2,... 

In this case it is also easy to determine the probabilities p n . We split the Poisson stream 
with rate A into the countably many independent Poisson streams numbered 1,2,.... The 
intensity of stream % is Ag« and the customers of this stream have a constant service time 
bi. The number of servers is infinite and thus we immediately obtain from the previous 
example that the number of type i customers in the system is Poisson distributed with 
mean Xqibi and this number is of course independent of the other customers in the system. 
Since the sum of independent Poisson random variables is again Poisson (see exercise 11) it 
follows that the total number of customers in the system is Poisson distributed with mean 
Xqibi + Xq 2 b 2 -\ = p. 

The examples above suggest that also in the M/G/oo the number of customers in 
the system is Poisson distributed with p. In fact, since each distribution function can be 
approximated arbitrary close by a discrete distribution function, this follows from example 
11.1.3 (see also exercise 74). Summarizing we may conclude that in the M/G/oo it holds 
that 

p n 
Pn = —re p , n = 0,1,2,..., 

n\ 

where p = XE(B). Note that this is true regardless of the form of the distribution function 

Fb(-) of the service time. 
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11.2 M/G/c/c queue 



In this model customers also arrive according to a Poisson process with rate A. Their 
service times are independent and identically distributed with some general distribution 
function. There are c servers available. Each newly arriving customer immediately goes 
into service if there is a server available, and that customer is lost if all servers are occupied. 
This system is therefore also referred to as the M/G/c loss system. 

In this section we want to find the probabilities p n of n customers in the system. Of 
special interest is the probability p c , which, according to the PASTA property, describes 
the fraction of customers that are lost. 

Example 11.2.1 (M/M/c/c) 

In this model the service times are exponentially distributed with mean 1/jj. 



A 
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Figure 11.2: Flow diagram for the M/M/c/c model 

From figure 11.1 we immediately obtain that 
Pn-iA = PnW, n = l,2,...,c. 
Hence, it is readily verified that 

(A//i) n /n! p n /n\ 



Pn J2 c n=0 (X/n) n /n\ j: c n=0 p n /n\ 
where p = A///. 



n = 0,l,...,c, 



It can be proved that also for a general service time distribution the probabilities p r< 
are given by (see e.g. [ ]) 



Pn 



P n l 



in 



71= 0,1,..., C, 



Z c n=oP n /n\> 

where p = XE(B). Hence, the so-called blocking probability B(c,p) is given by 

P c /c\ 



B(c,p) =p c 



2_/„=o P n I 



n\ 



11.2) 



The name given to this formula is Erlang's loss formula. Note that by Little's law we 
obtain that 

E(L) = p(l-B(c,p)) 
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Example 11.2.2 (Parking lot) 

Customers arrive according to a Poisson process at a parking lot near a small shopping 
center with a rate of 60 cars per hour. The mean parking time is 2.5 hours and the parking 
lot offers place to 150 cars. When the parking lot is full, an arriving customer has to park 
his car somewhere else. Now we want to know the fraction of customers finding all places 
occupied on arrival. 

The parking lot can be described by a M/G/150/150 model with p = 60 • 2.5 = 150. 
Hence the fraction of customers finding all places occupied on arrival is given by 

150 150 /150! 

B 150, 150 = ,. n ' — . 

1 ' £^150"/^ 

It will be clear that the computation of 5(150, 150) gives rise to a serious problem! 

11.3 Stable recursion for B(c,p) 

Fortunately it is easy to derive a simple and numerically stable recursion for the blocking 
probabilities B(c,p). From (11.2) we have 

B(c,p) = ^J^ 



Y, c n = P n ln\+p c /c\ 
Dividing the numerator and denominator of this expression by X^=o p n /n\ yields 

B(c,p) = ^V 1 '^ = ggfe-ZM . ( ii.3) 

K,HJ l + pB(c-l,p)/c c + pB(c-l,p) l ; 

Starting with B(0,p) = 1 we can use relation (11.3) to subsequently compute the blocking 
probabilities B(c, p) for c = 1, 2, 3, . . .. 



c 


£(c,150) 


150 


0.062 


155 


0.044 


160 


0.028 


165 


0.017 


170 


0.009 



Table 11.1: Blocking probability B(c,p) for p = 150 and several values of c 

Example 11.3.1 (Parking lot) 

Using the recursion (11.3) it is easy to compute B(c, 150) for c = 150. In table 11.1 we 
list B(c, 150) for several values of c. Clearly in the present situation with 150 parking 
places 6% of the arriving customers find all places occupied. With 20 additional places 
this percentage drops below 1%. 
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Remark 11.3.2 (Delay probability in the M/M/c queue) 

It will be clear from (5.1) that the formula for delay probability Hw in the M/M/c suffers 
from the same numerical problems as (11.2). Luckely there is a simple relation between 
the queueing probability HV in the M/M/c with server utilization p and the blocking 
probability B(c,cp). From (5.1) we immediately obtain 



n 



w - j--. * „ c _i 



(cp) c /c\ pB(c-l,cp) 



(1 - p) J2 c n -J (cp) n /nl + (cpY/c\ l-p + pB(c-l,cp) 

By first computing B(c — l,cp) from the recursion (11.3) we can use the relation above to 
determine U w . 

11.4 Java applet 

There is a JAVA applet is avalaible for the evaluation of the M/G/oo queue. The WWW- 
link to this applet is http://www.win.tue.nl/cow/Q2. 
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11.5 Exercises 

Exercise 74. 

Prove that the fact that (11.1) holds for discrete service time distributions implies that it 

also holds for general service time distributions. 

(Hint: Approximate the service time distribution from below and from above by discrete 

distributions and then consider the probability that there are n or more customers in the 

system.) 

Exercise 75. 

In a small restaurant there arrive according to a Poisson proces on average 5 groups of 

customers. Each group can be accommodated at one table and stays for an Erlang-2 

distributed time with a mean of 36 minutes. Arriving groups who find all tables occupied 

leave immediately. 

How many tables are required such that at most 7% of the arriving groups is lost? 

Exercise 76. 

For a certain type of article there is a stock of at most 5 articles to satisfy customer 
demand directly from the shelf. Customers arrive according to a Poisson process at a rate 
of 2 customers per week. Each customer demands 1 article. The ordering policy is as 
follows. Each time an article is sold to a customer, an order for 1 article is immediately 
placed at the supplier. The lead time (the time that elapses from the moment the order 
is placed until the order arrives) is exponentially distributed with a mean of 1 week. If on 
arrival of a customer the shelf is empty, the customer demand will be lost. The inventory 
costs are 20 guilders per article per week. Each time an article is sold, this yields a reward 
of 100 guilders. 

(i) Calculate the probability distribution of the number of outstanding orders. 

(ii) Determine the mean number of articles on stock. 

(iii) What is the average profit (reward - inventory costs) per week? 

Exercise 77. 

In our library there are 4 VUBIS terminals. These terminals can be used to obtain infor- 
mation about the available literature. If all terminals are occupied when someone wants 
information, then that person will not wait but leave immediately (to look for the required 
information somewhere else). A user session on a VUBIS terminal takes on average 2.5 
minutes. Since the number of potential users is large, it is reasonable to assume that users 
arrive according to a Poisson stream. On average 72 users arrive per hour. 

(i) Determine the probability that i terminals are occupied, i = 0, 1, . . . , 4. 

(ii) What is the fraction of arriving users finding all terminals occupied? 

(iii) How many VUBIS terminals are required such that at most 5% of the arriving users 
find all terminals occupied? 

116 



Exercise 78. 

Consider a machine continuously processing parts (there is always raw material available). 
The processing time of a part is exponentially distributed with a mean of 20 seconds. A 
finished part is transported immediately to an assembly cell by an automatic conveyor 
system. The transportation time is exactly 3 minutes. 

(i) Determine the mean and variance of the number of parts on the conveyor. 

To prevent that too many parts are simultaneously on the conveyor one decides to stop 
the machine as soon as there are N parts on the conveyor. The machine is turned on again 
as soon as this number is less than N. 

(ii) Determine the throughput of the machine as a function of N. 

(iii) Determine the smallest N for which the throughput is at least 100 parts per hour. 

Exercise 79. 

A small company renting cars has 6 cars available. The costs (depreciation, insurance, 
maintenance, etc.) are 60 guilders per car per day. Customers arrive according to a Poisson 
process with a rate of 5 customers per day. A customer rents a car for an exponential time 
with a mean of 1.5 days. Renting a car costs 110 guilders per day. Arriving customers for 
which no car is available are lost (they will go to another company). 

(i) Determine the fraction of arriving customers for which no car is available. 

(ii) Determine the mean profit per day. 

The company is considering to buy extra cars. 

(iii) How many cars should be bought to maximize the mean profit per day? 
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Solutions to Exercises 

Exercise 1. 

(i) Use that 



P(Y n > x) = n P(X { > x) 

1=1 

and 

n 

P(Z n <x) = l[P(X l <x). 

1=1 

(ii) It follows that 

P(X i = mm(X 1 ,...,X n )) = / "[[P(X j >x)fx l (x)d. 



t i i e-^ + -+^> dx = , ^ r. 

(yUl H V fJ-n) 



Exercise 1 
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Exercise 2. 

Use Laplace-Stieltjes transforms to prove that 

oo 

E(e~ sS ) = J2 P ( N = k ) E ( e ~ sS \ N = k ) 



k=i 

oo 



E(i-p)p a 



k=l 



K H + sJ n(l - p) + s 
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Exercise 2 



Exercise 3. 

Use that the random variable 

1/yUi, with probability pi, 
X = \ \ \ 

1/fik, with probability pk, 

has variance > 0, and hence that 

k / k \ ^ 



Exercise 3 
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Exercise 4. 

(i) Po (t) = P(A 1 >t) = e~ xt . 
(ii) Use that 

p n (t + At) = \Atpn-xit) + (1 - XAt)p n (t) + o(At), 

and let At tend to zero. 
(iii) Prove by induction that 

Pn(t) = K —^e~ xt 
n\ 

is the solution of the differential equation in (ii). 

Exercise 4 
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Exercise 5. 

(i) Po (t) = P(A 1 >t) = e~ xt . 
(ii) Use that 



P(N(t) =n)= P{N{t) =n\A 1 = x)f Al(x) dx. 
Jo 

(iii) Prove by induction that 

/ n ( Xt ) H At 

Pn(t) = ^e~ xt 
is the solution of the integral equations in (ii) . 



Exercise 5 
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Exercise 6. 

Merging property: Use that the minimum of independent exponential random variables is 
again an exponential random variable. (See also Exercise 1) 
Splitting property: Use the result of Exercise 2. 

Exercise 6 
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Exercise 7. 

Choose p = (18 - 4\/l4)/25 = 0.1213 and // = (2 - p)/4 = 0.4697. 

Exercise 7 
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Exercise 8. 

(i) For a Coxian-2 distribution we have 

£(X 2 ) = 4, + ^ + %, and E{X) = ± + *. 

Use this to show that E(X 2 ) > \E(X) 2 and hence that c\>\. 

(ii) Show that both distributions have the same Laplace-Stieltjes transform. Try to 
understand why these distributions are equivalent! (cf. Exercise 9) 

Exercise 8 
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Exercise 9. Use Laplace- Stieltjes transforms, or use the formula 

X x = mm{X x ,X 2 ) + (X, - mm{X 1 ,X 2 )) , 

where X\ is an exponential random variable with parameter A and X 2 is an exponential 
random variable, independent of Xi, with parameter /j — A. Exercise 9 
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Exercise 10. Use Exercise 9 with \i = \x\ and A = \x-z- Exercise 10 
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Exercise 11. Use generating functions and the fact that for the sum Z = X + Y of 

independent discrete random variables X and Y, it holds that (see subsection 2.2) 

Pz(z) = P x (z)-P Y (z). 

Exercise 11 
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Exercise 12. As time unit we take 1 minute. 
(i) Solve the (global) balance equations 

MnPn = APn+1, 71=0,1,2,3, 

where A = \i = 1/3, together with the normalization equation. This gives 

32 32 24 12 3 

P0= 103' Pl= 103' P2= 103' P3= 103' P4 = To3- 

(ii) E(L) = 128/103 « 1.24. 

(hi) £7(5) = 384/71 » 5.41 minutes. 

(iv) £7(5) = 384/103 « 3.73 minutes. 
£(W) = 171/103 » 1.66 minutes. 

Exercise 12 
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Exercise 13. 

(i) Exponential with parameter p* = p(l — p) (see Exercise 2). 
(ii) P(L = n) = (1 - p)p n for n = 0, 1, 2, ... , where p = A///*. 

Exercise 13 
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Exercise 14. We have that 
service completion rate 



Hi, if nr. of customers < Ql, 

/x, if Ql < nr. of customers < Qh, 

Hh, if nr. of customers > Qh- 



The (global) balance equations are 

Ap n = VLPn+i, ifn + l<Q L , 

Xp n = np n+1 , if Q L < n + 1 < Q H , 

Apn = HHPn+i, ifn + l>Q H - 

The solution of these equations is given by 



Pn 



A^-^A 



P() 

Po 
Finally, po follows from the normalization equation. 



,+i 



-Qi 



k-Qh+1 



if n < Ql, 

if Ql < n < Qh, 

if n > Q H - 



Exercise 14 
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Exercise 15. 

(i) 3/8 
(ii) 5/3 
(iii) 80 minutes 

Exercise 15 
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Exercise 16. 

(i) p( L = n ) = l(!) n 7 n = 0,l,2,... 

and hence E(L) = 2 and cr 2 (L) = 6. 

(ii) P(S < t) = 1 - e" 4 / 6 , £ > 0, 
P(W < t) = 1 - fe"*/ 6 , t > 0. 

(iii) fe" 1 / 3 « 0.48. 

(iv) Po = §, Pi = §, P2 = j|, 

hence £(L) = 14/19 » 0.737 and a 2 (L) = 22/19 - (14/19) 2 » 0.615. 

(v) £(£) = 42/19 « 2.21 minutes and £(W) = 12/19 » 0.63 minutes. 

Exercise 16 
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Exercise 17. 

(i) It holds that 



and 



1 /9\ n 

P(L«*») = n ) = - (_j } n = ,l,2,... 



P(L (LPG) = n) = 5 -Q n , n = 0,l,2,.... 



(ii) Use 



P(L = n ) = £ p(l,( G <") = fc, L( LPG ) =n-k) 

k=0 



to show that 



p( L = n) = — (-)"- — (-) . n=0.L2.... 
V ; 54 V3/ 54 V6/ 
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Exercise 17 



Exercise 18. 




(i) E(S 1 ) = 


■ 7.5 minutes. 


E(S 2 ) = 


: 30 minutes. 


(ii) E(St) = 


: 10.625 minutes 


E(S 2 ) = 


: 27.5 minutes. 



Exercise lc 
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Exercise 19. Fraction of time that it is crowded: p 5 ~ 0.24. 

Number of crowded periods (per 8 hours = 480 minutes): Xp^ = 480(1 — p)p 4 ~ 38. 

E (crowded period) = E(busy period) = 3 minutes. Exercise 19 
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Exercise 20. We have 

average costs per hour = 16/z + 20E(L q ) 

P 2 
= 16// + 20 ' 



1-P) 

8000 
16/t 



//(// - 20) 
For // > 20, this function is minimal for \x = /** « 25. Exercise 20 
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Exercise 21. As state description we use the number of jobs in the system, and if this 
number of jobs is equal to 1, we distinguish between state (1, /) in which the fast server is 
working and state (1, s) in which the slow server is working. 

(i) As solution of the balance equations we find 

Pl = PlJ + Pl,s = Cpo , 
Pn = P u ~ 1 Pi , n> 1, 

where we used the notation /i = \x\ + p,2, P = A/// and 

A^(A + /i 2 ) 
/ii/i 2 (2A + /i)' 

(ii) For the mean number of jobs in the system we find 

C 



E { L ) = J2 n Pn = J 1 



n=l 



(l-p)(l-p + C) 



(iii) It is better not to use the slower machine at all if E(L^), the expected number of 
jobs in the system when you only use the fast server is smaller than E(L). This is 
the case if \i\ > A and 



A C 

< 



(jh-x - (i-p)(i- P + cy 

(iv) In case (a) we have 

E(Lf) = - < — - E(L). 
1 ; 3 104 v ; 

In case (b) we have 

E(L') = \>f 7 = E{L). 



Exercise 21 
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Exercise 22. As time unit we choose 1 minute: A = 4/3 and \x = 1. In order to have 
p < 1 we need that c > 2. Hence, we first try c = 2. This gives Uw = 8/15 ~ 0.533 and 
E(W) = 24/30 = 0.8 minutes. Hence, we conclude that 2 boxes is enough. Exercise 22 
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Exercise 23. As time unit we choose 1 minute: A = 2/3 and // = 1/3. In order to have 
p < 1 we need that c > 3. Hence, we first try c = 3. This gives Tiw = 4/9 ~ 0.444 and 

P(W > 2) = Il w ■ e" 2/3 w 0.228 > 0.05. 
Similarly, for c = 4 we find II ^ = 4/23 w 0.174 and 

P(H/ > 2) = n w • e" 4/3 w 0.046 < 0.05. 
Hence, we need at least 4 operators. Exercise 23 
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Exercise 24. As time unit we choose 1 minute: A = 1/3 and p = 1/3. 
(i) We have 

1 1 /l\ n_1 

Po = ~ s , Pn = - s {- 2 ) , n>l. 

(ii) Using (5.2) and (5.3) we have 

E(L q ) = U w ■ -^- = 1/3, E(W) = U w ■ - 1 — ■ — = 1 minute. 
\ — p 1 — p Cfl 

(hi) U w = 1/3. 

(iv) For c = 2 we have TL\y = 1/3 > 1/10. Similarly, we can hnd for c = 3 that 
Tlw = 1/11 < 1/10. Hence, we need 3 troughs. 

Exercise 24 
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Exercise 25. As time unit we choose 1 minute: A = 15 and \x = 6. 
(i) c - p = A/// = 2.5. 

(ii) c • (1 — p) ■ 12 = 6 maintenance jobs per minute, 
(iii) We have 

8 20 25 /5\ n " 2 



P0= 178' W= 178' P - = 178 U; ' n " 2 ' 

and hence (see (5.1)) 

125 

n w = « 0.702. 

w 178 

(iv) Using (5.3), we have 

£(W) = n w = - • Tlw ~ 0.234 minutes. 

1 — p cfi 3 



Exercise 25 
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Exercise 30. 

(i) The distribution of the number of uncompleted tasks in the system is given by 

7 /2\ n 7 ( 2\ n 
""=24(3) + io(-5) ' » = <U*.-- 

(ii) The distribution of the number of jobs in the system is given by 

35 /4\ n 21 / 4 \ n 
fe= 48UJ "SOU; ' ^ = 0,1,2,.... 

(iii) The mean number of jobs equals E(L) = J2^Li n Qn = 104/105. 
(iv) The mean waiting time of a job equals equals 

E(W) = -^—E(R B ) = 8/7 • 3/2 = 12/7 minutes. 
1-p 

(Check: Little's formula E(L) = XE(S) is satisfied, with E(L) = 104/105 job, A = 4/15 
job per minute and E(S) = 26/7 minutes.) Exercise 30 
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Exercise 31. Define Tj as the mean time till the first customer is rejected if we start with 
i phases work in the system at time t = 0. Then we have 



To = 


= 1 + T 2 , 


Ti = 


1 1 1 

= 2 + 2 To+ 2 Ts ' 


T 2 = 


1 1 1 
" 2 + 2 T ' + 2^ 


T 3 = 


1 1 

" 2 + 2 T - 


T 4 = 


" 5 + 5 T - 



The solution of this set of equations is given by 
(T ,T x ,T 2 ,T 3 ,Tt) = {4:,3±,3,2,l±). 

Hence, if at time t = the system is empty, the mean time till the first customer is rejected 
is equal to 4. Exercise 31 
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Exercise 32. The distribution of the number of phases work in the system is given by 

7 (2\ n 7 ( 2\ n 
""=24(3) + io(-5) ' » = <U.*-- 

(i) The distribution of the waiting time (in minutes) is given by 

7 1 1 7 

v ~ ' 12 20 

(ii) The fraction of customers that has to wait longer than 5 minutes is given by 

7 5 1 7 

P(W > 5) = —e~^ e" « 0.376. 

v ; 12 20 

Exercise 32 



150 



Exercise 33. 

(i) The distribution of the number of customers in the system is given by 

3 /l\ n 6 / l\ n 
""=7(2) + 35 Hi) ' » = M,2,.... 

(ii) The mean number of customers equals E(L) = J2^=i np n = 5/6. Now, either use the 
PASTA property 

E(S) = (5/6) -6 + (1/4) -6 + 6 

or use Little's formula 

5/6 



E(S) 



1/15 



to conclude that the mean sojourn time of an arbitrary customer is equal to 12.5 
minutes. 

Exercise 33 
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Exercise 34. 

(i) See Example 6.2.1. 

(ii) Use PASTA and/or Little to conclude that E(S) = 11/12 week. 

(iii) Because po + Pi + Vi < 0-99 and Po + Pi + P2 + P3 > 0.99 we need at least 4 spare 
engines. 

Exercise 34 
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Exercise 35. 

(i) The distribution of the number of uncompleted tasks at the machine is given by 

(ii) The mean number of uncompleted tasks equals E(L task ) = J2™=i n Pn = 11/4. Hence, 
using PASTA, the mean waiting time of a job is E(W) = E(L task ) ■ 1 = 11/4 minutes. 

(iii) The mean sojourn time of a job equals equals E(S) = E(W) + E(B) = 11/4 + 8/5 = 
87/20 minutes. Hence, using Little's formula, we have E(Lj b) = 5/12-87/20 = 29/16. 

Exercise 35 
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Exercise 36. 

(i) The distribution of the number of customers in the system is given by 

p - = 5 (2) + T5h) • » = m.*.-- 

(ii) For the mean number of customers in the system we have E(L) = S^Li nPn = 3/4. 
Hence, using PASTA, the mean waiting time of the first customer in a group equals 
E(Wi) = E(L) ■ 5 = 15/4 minutes. 

(iii) The mean waiting time of the second customer in a group equals E{Wz) = E(W\) + 
5 = 35/4 minutes. 

(Check: Little's formula E(L q ) = XE{W) is satisfied, with E(L q ) = 3/4 - 1/3 = 5/12 
customer, A = 1/15 customer per minute and E(W) = 25/4 minutes.) Exercise 36 
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Exercise 37. As time unit we take 1 minute. Hence, A = 1/2 and 

1 i 3 1 



4 i+s 4 1 + s 



B(s) 

(i) From (7.6) we have 

(l-p)B(A- Xz)(l-z) 3 15-72 



Pl(z) 



9_ 
32 



3_ 
32 



B(X-Xz)-z 8 (3-2z)(5-2z) 

(ii) The distribution of the number of customers is given by 

9 /2\ n 3 {2\ n 

, n = 0,1,2,.... 



Pn 32 V37 ' 32 V5 
(hi) £(L) = £- =1 np n = 43/24 
(iv) From (7.7) we have 



S(s) 



(1 - p)B(s)s S 4 + 7s _ 27 

XB(s) + s - A ~ 4 (1 + 4s)(3 + 4s) ~ 32 



(v) The distribution function of the sojourn time is given by 



F s (x) 



27 



(l-e"i s 



5 



32 v '32 

and the mean sojourn time by 

43 



(1 



e 4 



'), 



y J 32 32 3 

(vi) From (7.21) we have 



12 



minutes. 



, s E(B) 10 
E(BP) = -^-^- = — minutes. 
1 — p 3 

(vii) For the M/M/l queue we have 



and 



E(L) 



E(S) 



5 



1-p 3' 



E(L) _ 10 
A ~ T 



minutes. 



5 
32 



Exercise 37 
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Exercise 38. As time unit we take 1 minute, so A = 1/6. 



CO B(a) = (+ 



2 



(ii)Pn = f(|) n -^(|) n , n = 0,l,2,.... 

(iii) £(L) = 11/24 and E(S) = 11/4 minutes. 

Exercise 38 
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Exercise 39. As time unit we take 1 minute, so A = 1. Let X be exponential with 
parameter 4 and Y be exponential with parameter 1. Then, 

11 1 3 

E(B) = - ■ E(- + X) + -- E(Y) = - minutes, 

v i 2 4 2 4 

E(B<) - I. B([ i + ^) + i-^) 444-2= | minutes, 

and so E(R) = 37/48 minutes. Hence, E(W) = 37/16 minutes and E(L q ) = 37/16 
customers. Exercise 39 
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Exercise 40. As time unit we take 1 minute, so A = 1/20. Furthermore, E(B) = 12 
minutes and E(R) = 31/3 minutes. Hence, E(S) = 55/2 = 27.5 minutes. Exercise 40 
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Exercise 41. The mean waiting time of jobs is given by 

o 25 

E(W) = ^— ■ E(R) = — minutes. 
1 - p 4 

Exercise 41 
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Exercise 44. The time unit is 1 minute: A = 1/6, E(B)=15/A, p = 5/8, E(R) = 
(2/5) -6 + (3/5) -3 = 21/5. 

(i) The service time is hyperexponentially distributed with parameters p\ = 1/4, p2 = 
3/4, fii = 1/6 and /x 2 = 1/3. 

(ii) From (7.9) we have 

~, , (1-P)s l + 9s + 18s 2 3 9 1 11 

Wis) = —=> — = = - H 1 . 

XB(s) + s-X (l + 12s)(l + 4s) 8 16 1 + 12s 16 1 + 4s 

(iii) The distribution function of the waiting time is given by 

Fw(x) = | + A (1 - e -*«) + 1 (1 - e"h- 
Hence, the fraction of cows for which the waiting time is less than 3 minutes equals 
Fw(S) = | + A (1 - e-*) + ^ (1 - e-f) « 0.532. 

(iv) The mean waiting time is given by 

9 1 

E(W) = — ■ 12 H 4 = 7 minutes. 

16 16 

Alternatively, from a mean value analysis we have 

£7(W) = —^E(R) = = 7 minutes. 

I — p 3 5 

Exercise 44 
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Exercise 45. The time unit is 1 hour: A = 1, E(B)=7/12, p = 7/12, E(R) = (3/7) 
(7/12) + (4/7) • (1/3) = 37/84. 

(i) The Laplace-Stieltjes transform of the processing time is given by 

s, n 4 3 

B(s) 



4+s 3+s 
From (7.7) we have 

(l-p)B(s)s 



S(s) 



XB(s) + s-X (l + s)(5 + s) 4 1 + s 4 5 + s 
(ii) The distribution function of the production lead time is given by 

F 5 (x) = ^(l-e-)-i(l-c- fa ). 

The mean production lead time is given by 

E(S) = --1---- = - hours. 
v ' 4 4 5 5 

Alternatively, from a mean value analysis we have 

E(S) = r 4-E( fi ) + B ( B) = Z.g + i U5 hours. 



(iii) The mean cost per hour equals 
A • (1 - F s (3)) ■ 100 = ( 



5 1 

- ■ e -3 e -15 ) • 100 « 6.22 dollar. 

4 4 



Exercise 45 
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Exercise 46. The time unit is 1 minute: A = 1/10, E(B)=25/4, p = 5/8, E(R) = 
(2/5) -10 + (3/5) -5 = 7. 

(i) The pick time is hyperexponentially distributed with parameters p\ = 1/4, p2 = 3/4, 
Hi = 1/10 and \i2 = 1/5. 

(ii) From (7.7) we have 

(l-p)B{s)s 12 + 105s 5 3 27 1 



S(s) 



\B(s) + s-X 4(3 + 20s)(l + 20s) 32 3 + 20s 32 1 + 20s 
(iii) The distribution function of the sojourn time is given by 

^(x) = |(l-e-^) + |(l-e-^). 

Hence, the fraction of orders for which the lead time is longer than half an hour is 
given by 

, s 5 9 27 3 

1 - F s (30) = — ■ e~i + — • e"* w 0.190. 
v ; 32 32 

(iv) The mean lead time is given by 

_, m 5 20 27 nn 215 . 

E(S) = 1 20 = minutes. 

v ; 32 3 32 12 

Alternatively, from a mean value analysis we have 

E(S) = ^-E(R) + E(B) = - ■ 7 + - = — minutes. 
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Exercise 51. As time unit we choose 1 minute: \x = 1. The Laplace-Stieltjes transform 
of the interarrival time distribution is given by 

A(s) 



3 1 + s 3 l + 3s 
(i) a n = (1 — a) a n , where a, the solution in (0,1) of a = A(n — fj,a), is given by 

0.48. 





21 


- V153 






18 


(ii) E(L a ) - 


_ O" 
1-CT 


« 0.92. 


(hi) S(s) = 


l-O" 
1-CT+S' 




(iv) E(S) = 


^- » 1.92. 

1— <T 



(v) £(L) = A • E(S) = f • E(S) « 0.82 
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Exercise 52. We have yu = 6 and the Laplace-Stieltjes transform of the interarrival time 
distribution is given by 

~, . 13 3 11 2 
A(s) = — 1 . 

v ' 24 3 + s 24 2 + s 

(i) a n = (1 — a) d", where a, the solution in (0,1) of a = A(/j — /j<j), is given by a — '' 



12" 



(ii) F w (t) = 1 - a e-^ 1 -^ = l _ A e -| * 



12 
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Exercise 53. The sojourn time is exponentially distributed with parameter fi(l — a), 
where \x = 1 and 



Exercise 53 



7- V17 
a = ^ — w 0.36. 
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Exercise 54. 

(i) The solution in (0,1) of a = e~ 2 ( 1_cr \ is given by a « 0.203. 
(ii) F s (t) = 1 - e-^ 1 -^ « 1 - e -(°- 4 )-*. 
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Exercise 55. 

(i) 3/8. 
(ii) Let p n denote the probability that there are n cars waiting for the ferry. Then, 

1 3 /1\ 2 7 



Po = — — , 

y 4 4 \2J 16' 

3 /IN 1 3 /1\ 3 15 



Pl 4 V27 '4 V2/ 32' 

3 /l\ n+2 
Pn = r y - ">2- 

(iii) E(L q ) = Yln°=o n Pn = f ? an d hence using Little's formula we have -E^W) = 3 minutes. 
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Exercise 56. As time unit we choose 1 minute: 
A = 1/6, E(B) = 9/2, p = 3/4 and E(R) = 5, 

Ai = 1/12, E(Bi) = 3,pi = 1/4 and E(R l ) = 3, 

A 2 = 1/12, E{B 2 ) = 6,p 2 = 1/2 and £(i? 2 ) = 6. 

(i) E(W) = j^- p E{R) = 15 minutes, 
(ii) Use formula (9.3) on page 89: 

p 1 E(R 1 ) + p 2 E(R 2 ) . 

-E/fW-i) = = 5 minutes , 

1 ~ Pi 

Pl E(R 1 ) + p 2 E(R 2 ) . 

#(Vv 2 ) = 7j rr; r = 20 minutes , 

(i-Pi)(i-Pi-p2j 



E(W) = 1 E{W 1 ) + -E(W 2 ) = 12.5 minutes . 

(iii) Similar to formula (9.3), we now have 

Pl E{R{) + p 2 E{R 2 ) 15 . 

£/(vr 2 ) = = — = 7.5 minutes , 

1 - p 2 2 



Pl E{R{) + p 2 E{Ih) „ n . . 

£ (Wi) = 7- ry- r = 30 minutes , 

(l-p 2 )(l-pi - p 2 ) 



E(W) = -£(Wi) + ]rE(W 2 ) = 18.75 minutes 
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Exercise 57. As time unit we choose 1 minute: 

Ai = 1/60, £(£i) = 10, pi = 1/6 and £(#i) = 5, 

A 2 = 1/30, E(B 2 ) = 10, p 2 = 1/3 and £(i? 2 ) = 5, 

A 3 = 1/30, E(B 3 ) = 10, & = 1/3 and E(R 3 ) = 5. 

Now, use formula (9.5) for E(Si) on page 90: 

= gigCgO + E , B j = n minutes 

p 1 E(R 1 ) + p 2 E(R 2 ) E(B 2 ) . 

£5 2 = — + = 18 minutes , 

(1-Pl)(l-Pl -P2j 1-Pl 

P,, q x PlEjRJ + p 2 E(R 2 ) + p 3 E(R 3 ) gCg3J yn . . 

E(S 3 ) = — + = 70 minutes , 

(1 -pi -p 2 )(l - Pi - P2 - p3) I-P1-P2 
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Exercise 58. As time unit we choose 1 minute: 

A = 1/15, E{B) = 55/4, p = 11/12 and E{R) = 15/2, 

Ai = 1/30, £(£i) = 10, pi = 1/3 and E(R l ) = 5, 

A 2 = 1/60, £(£ 2 ) = 15, p 2 = 1/4 and E(R 2 ) = 15/2, 

A 3 = 1/60, -E(£ 3 ) = 20, p 3 = 1/3 and E{R 3 ) = 10. 

(i) E{W 1 ) = E(W 2 ) = E(W 3 ) = E(W) = j^E(R) = 165/2 = 82.5 minutes. Hence, 
E(Si) = 92.5 minutes, E(S 2 ) = 97.5 minutes, E(S 3 ) = 102.5 minutes and E(S) = 
96.25 minutes. 

(ii) Use formula (9.4) for E(Si) on page 89: 

_ Pl E(R l) + P2 E(R 2 ) + P3 E(R 3 ) + _ 325 ^ M ^^ 

1-pi 16 

p ,,, p^jR^ + p 2 E(R 2 ) + p 3 E(R 3 ) , „ m . 159 . 

E(5 2 ) = 7- rr- : h -E(-B 2 ) = — — = 39.75 minutes , 

(l-pi)(l-pi-p 2 ) 4 

mq\ PiE{R{) + p 2 E{R 2 ) + p 3 E(R 3 ) ^^ . 

E(5 3 ) = 7- — r + E(B 3 ) = 218 minutes , 

(l - pi - p 2 )(l - pi - p 2 - p 3 ) 

E(S) = )rE{S{) + jE{S 2 ) + jE{S 3 ) « 74.59 minutes . 



(iii) Combine the arguments of Sections 9.1 and 9.2: 



245 

E(Si) = — ^r ' ^ V ^ Z ' + E(B 1 ) = « 15.31 minutes 

1 16 



£(S 3 



Pl£(fll) + 


p 2 ^(i? 2 ) 


-E(B t ) = 


1 - 


Pi 


Pl£(i?l) 4 


- p 2 £(i? 2 ) 


+ £(£ 2 ) : 


(1-Pl)(l 


- Pi - P2) 


PiE{R{) + 


p 2 E{R 2 ) + p 3 £(i? 3 ) 



£(£2) = /,'" '"!,',' r ^ "•■■■ '-. +_ £(# 2 ) = = 27.75 minutes , 



E(B 3 



;i -pi -p 2 )(i -pi -p 2 -p 3 ) 1-P1-P2 



246 minute s , 



E(S) = -E(Sx) + -£(<%>) + -E(S 3 ) » 76.09 minutes 
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Exercise 59. As time unit we choose 1 minute: A = 1/6. 

(i) For N, the number of parts that has to be produced for an order, we have 

P(N = n)= Q\ n=l,2,3,... 

and hence E(N) = 2 and a 2 (N) = 2 (see also Section 2.4.1). From B = 2N, it now 
follows that E(B) = 4 and a 2 (B) = 8. 

(ii) Using that p = 2/3 and E(R) = 3, we have 

E(S) = -?—E(R) + E(B) = 10 minutes . 
1-p 

(iii) We now have 

Ai = 1/12, ^(J3i) =2,pi = 1/6 and E{R 1 ) = 1, 

A 2 = 1/12, E{B 2 ) = 6,p 2 = 1/2 and E{R 2 ) = 11/3. 

Hence, 

_ MlftltMW + E(Bi) _ 22 _ 4 4 m . nutes 

1 -Pi 5 

PiEjRj + p 2 E(R 2 ) 66 . 

^(5 2 ) = t- r-. r + E(B 2 ) = — = 13.2 minutes , 

(1-Pi)(l-Pi-P2j 5 

(iv) E(S) = \E(Si) + \E(S 2 ) = f = 8.8 minutes. 

Exercise 59 



171 



Exercise 63. 



E ( S) = ^E(R B ) + ¥ ^r T - + E(B). 
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Exercise 64. 

As time unit we choose 1 second: 
A = 1/6. 

(i) The mean waiting time satisfies 

E(W) = 2.5 + £(L 9 ) • 5 . 

Together with Little's formula, E(L q ) = XE(W), this yields E(W) = 15 seconds. 

(ii) The time elapsing from entering the carrier till the departure of that bin is 4 cycles 
(= 4 • 5 = 20 seconds) plus moving out of the carrier (= 2 seconds), so 22 seconds. 
Hence, the mean sojourn time is equal to 15 + 22 = 37 seconds. 
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Exercise 65. 

(i) The mean number of orders in the system is given by 

/ s P N-l 

E(L) = — — + . 

v ; 1-p 2 

(ii) From Little's formula we obtain 

1 — p 2A 
(iii) The average cost (setup cost + machine cost + waiting cost) per minute equals 



6 n ,. 3 /. r XX 6 21 3N 

— + 8 + 4 + -(N - 1)) = — + — + — . 
iV v 2 v ;; N 2 2 



(iv) N = 2. 
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Exercise 66. See exercise 4 of the exam of June 21, 1999. Exercise 66 
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Exercise 69. As time unit we choose 1 hour: 

A = 1,E(B) = 1/2, p = 1/2 and E(R B ) = 1/2. 

(i) The fraction of time that the machine processes orders is 1/2. The mean duration of 
a period that the machine is switched off equals 1, the mean duration of a switch-on 
period equals T. Hence, the mean duration of a period that the machine processes 
orders equals 1+T. Hence, both the mean number of orders processed in a production 
cycle and the mean duration of a production cycle equals 2 + 2T. The mean waiting 
time of an order equals 

1 T T 1 + T 

E(W) = E(L q ) ■ 1/2 + T + + • 1/2. 

v ' y ' ' 2 + 2T 2 + 2T 2 2 + 2T ' 

Together, with Little's formula E(L g ) = 1 • E(W) this gives 

E(w) = r 2 + 3r + i 

v ; 2 + 2T 

Hence, the mean production lead time of an order equals 

T 2 4- AT 4- 2 

E(S) = - 



2 4-2T 
(ii) The average cost per hour equals 

17 T 2 4- 3T 4- 1 T 2 + 3T + 18 



2 4- 2T 2 4- 2T 2 4- 2T 

which is minimal for T = 3. 
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Exercise 71. 

(i) E(S) = 105/2 = 52.5 minutes . 
(ii) E(L) = 21/6. 
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Exercise 73. As time unit we choose 1 minute: 

A = 1/10, E{B) = 15/2, p = 3/4 and E(R B ) = 35/9. 

(i) The fraction of time that the server serves customers is 3/4. The mean duration of 
a period that the server is away equals 10 + 10 + 5 = 25 minutes. Hence, the mean 
duration of a busy period equals 75 minutes. 

(ii) 75/7.5 = 10 customers. 

(iii) The mean waiting time of a customer equals 

E(W) = E{L q ) ■ 15/2 + 1/10 • 5 + 1/20 • 5/2 + 3/4 • 35/9. 

Together, with Little's formula E(L q ) = 1/10 • E(W) this gives E(W) = 121/6 = 
20.17 minutes. Hence, the mean sojourn time of a customer equals E(S) = 27.67 
minutes. 
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Exercise 77. 

(i) For the probability that i terminals are occupied we have 

!r 8 3 l 



Pi ^4 3^ 

■>n=Q n! 



Hence, 

. / 8 24 36 36 27 \ 
(P0,Pi,p 2 ,P3,P4) = (^ Wl , I g T , 131, Y^J • 

(ii) 5(4,3) =p 4 = ^ =0.2061. 
(hi) Use the recursion (11.3): 

5(4, 3) = 0.2061, 5(5, 3) = 0.11005, 5(6, 3) = 0.05215, 5(7, 3) = 0.0219. 

So, we need at least 7 terminals. 
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Exercise 79. 

(i) 5(6,7.5) = 0.3615. 
(ii) The mean profit per day equals 

5 • 110 • 1.5 • (1 - 5(6, 7.5)) - 6 • 60 = 166.7 guilders . 

(iii) When the company has c cars, the mean profit per day equals 

5 • 110 • 1.5 • (1 - B(c, 7.5)) - c • 60 guilders . 

So, if the company buys 1 extra car, the mean profit becomes 174.7 guilders, if the 
company buys 2 extra cars, it becomes 173.8 guilders, if the company buys 3 extra 
cars, it becomes 163.4 guilders, and so on. The mean profit per day is maximized 
when the company buys 1 extra car. 
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